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STABLE PATTERNS WITH JUMP DISCONTINUITY IN SYSTEMS 
WITH TURING INSTABILITY AND HYSTERESIS 


STEFFEN HARTING, ANNA MARCINIAK-CZOCHRA, AND IZUMI TAKAGI 

Abstract. Classical models of pattern formation are based on diffusion-driven insta¬ 
bility (DDI) of constant stationary solutions of reaction-diffusion equations, which leads 
to emergence of stable, regular Turing patterns formed around that equilibrium. In this 
paper we show that coupling reaction-diffusion equations with ordinary differential equa¬ 
tions (ODE) may lead to a novel pattern formation phenomenon in which DDI causes 
destabilization of both constant solutions and Turing patterns. Bistability and hysteresis 
effects in the null sets of model nonlinearities yield formation of far from the equilibrium 
patterns with jump discontinuity. We derive conditions for stability of stationary solu¬ 
tions with jump discontinuity in a suitable topology which allows us to include the dis¬ 
continuity points and leads to the definition of (eo, A)-stability. Additionally, we provide 
conditions on stability of patterns in a quasi-stationary model reduction. The analysis 
is illustrated on the example of three-component model of receptor-ligand binding. The 
proposed model provides an example of a mechanism of de novo formation of far from 
the equilibrium patterns in reaction-diffusion-ODE models involving co-existence of DDI 
and hysteresis. 


1. Introduction 

Since the seminal paper of Alan Turing [26], mathematical models of biological pat¬ 
tern formation have been constructed by using systems of reaction-diffusion equations 
exhibiting diffusion-driven instability (DDI). DDI is related to a local behavior of solu¬ 
tions of a reaction-diffusion system in the neighborhood of a constant stationary solution 
that is destabilized through diffusion. It may lead to emergence of stable continuous and 
spatially periodic structures around the destabilized constant equilibrium. The Turing 
concept became a paradigm for pattern formation and led to development of numerous 
theoretical models, though its biological verification has remained elusive [DEI El]. 

The models are based on the idea that cells differentiate according to positional in¬ 
formation which is supplied to them by diffusing biochemical morphogens. Different 
morphogen concentrations are able to activate transcription of distinct target genes and 
thus lead to cell differentiation. However, both regulatory and signaling molecules (lig¬ 
ands) act by binding and activating receptor molecules, which are located in the cell 
membrane HUES]. This observation leads to a hypothesis that the positional value of a 
cell may be determined by the density of bound receptors which do not diffuse. Taking 
dynamics of receptors into account leads to systems coupling reaction-diffusion equations 
with ordinary differential equations (ODE) [9] ED 13], called also receptor-based mod¬ 
els. Such systems can be obtained as a homogenization limit of the models describing 
coupling of cell-localized processes with ccll-to-cell communication via diffusion in a cell 
assembly mm- Nonlinear interactions of diffusive and non-diffusive components arise 
also from modeling of interactions between cellular or intracellular processes. Reaction- 
diffusion-ODE models have recently been employed in various biological contexts, see e.g., 

[D El ED ED ED ED- 
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Although receptor-based models may exhibit Turing instability as shown in ra and 
discussed more recently on several examples from mathematical biology in [9j, they are 
different from classical Turing-type models. Since DDI is induced by self-enhancement 
in the subsystem with smallest diffusion [2], in reaction-diffusion-ODE models instability 
is induced by growth properties of the ODE subsystem. Therefore, coupling of reaction- 
diffusion equations to ODEs may lead to DDI which does not arise if the non-diffusive 
components are neglected. 

To understand the role of non-diffusive components in pattern formation process, we 
focus on systems involving a single reaction-diffusion equation coupled to ODEs on one¬ 
dimensional domain I. In general, equations of such models can be represented by the 
following initial-boundary value problem, 


( 1 . 1 ) 


du 

dt 

dv 

dt 

(m(0),u(0)) 


d 2 v 

(uo,v 0 ) g C(J) X (C 2 (/)nC(7)), 


x e l, t > o, 
x & I, t > 0, 


with homogeneous Neumann boundary conditions for v. 

It is an interesting case, since a scalar reaction-diffusion equation cannot exhibit stable 
spatially heterogeneous patterns [3]. Coupling it to an ODE fulfilling autocatalysis condi¬ 
tion, i.e. df /du > 0 at the constant equilibrium, leads to DDI. However, in this case all 
Turing patterns are unstable, i.e. the same mechanism which destabilizes constant solu¬ 
tions, destabilizes also all continuous spatially heterogeneous stationary solutions [15, 16). 
This instability result holds also for patterns with jump discontinuity in case of a special 
class of nonlinearities [15], 16]. The question then arises as to which patterns, if any, can 
be exhibited in such models. 

As shown in na, it may happen that there exist no stable stationary patterns and the 
emerging spatially heterogeneous structures have a dynamical character. Simulations of 
different models of this form indicate formation of dynamical, multimodal and apparently 
irregular structures, the shape of which depends strongly on initial conditions m m m- 
On the other hand, reaction-diffusion-ODE models may give rise to discontinuous patterns 
due to the hysteresis effect in the model nonlinearities, i.e. when the equation f(u,v) = 0 
has multiple quasi-stationary solutions u; = Hi{u). Diffusion tries to average different 
states and may lead to formation of far from equilibrium structures. Hysteresis yields 
emergence of stationary solutions with a jump discontinuity, which may be monotone, 
periodic or irregular. The hysteresis-driven pattern formation has been investigated an¬ 
alytically in systems with two stable constant equilibria, i.e. without DDI [Qj [T5], TU] , 
In such a case patterns depend strongly on initial conditions and require large initial 
heterogeneity. 

Consequently, the next question is whether a coupling of the two mechanisms, DDI and 
hysteresis, can provide a mechanism of de novo formation of stable far from equilibrium 
patterns. By de novo pattern formation we understand emergence of spatially hetero¬ 
geneous structures starting from small perturbation of the constant stationary solution. 
According to our knowledge, systems with DDI and hysteresis have not been studied so 
far. 

Motivated by these observations, we propose a new receptor-based model exhibiting 
DDI and hysteresis which lead to de novo formation of stable patterns with jump-type 
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discontinuities. The model is a modification of the simplest model describing receptor- 
ligand binding dynamics developed in [19, TSj. The modification accounts for saturation 
in production terms, which provides uniform boundedness of the model solutions. We 
consider two- and three-component systems, i.e., 

• a model consisting of one ordinary differential equation and one reaction-diffusion 
equation (10DE+1RDE), 

• a model consisting of two ordinary differential equations and one reaction-diffusion 
equation (20DE+1RDE), 

While the principle of DDI is similar to Turing’s original idea, coupling to an ODE is 
crucial for stability of steady states in this paper. One interpretation of the observed 
phenomenon is that non-diffusive species, such as receptors, can facilitate formation of 
patterns with sharp transitions. 

In this paper, we derive conditions for stability of such patterns. Since stationary 
solutions for nondiffusing variables have jump discontinuities, linearized stability analysis 
requires careful treatment. To cope with this difficulty, we consider a special topology 
which allows us to include the discontinuity points and yields a notion of (e 0 , /^-stability. 
Our approach follows the one of Weinberger proposed in refs. [28, [3] to study discontinuous 
solutions in a model of density dependent diffusion. A novelty of our paper lies in the 
stability analysis of patterns with jump discontinuity in systems involving several ODEs 
which may be reduced to a subsystem by using quasi-steady state approximation. We 
show that under suitable conditions on the nonlinearities, stability of stationary solutions 
is preserved under quasi-steady state reduction. 

The paper is divided into five sections: In Section [2] the main results of this paper are 
formulated. We consider systems of type (11.111 with scalar-valued v and u having values in 
either M or K 1 2 and provide conditions under which steady states with jump discontinuity 
are stable in the topology introduced in ref. [28]. Then, we extend the stability results to 
systems with a small parameter and the corresponding quasi-stationary approximation. In 
Sectional we present a specific receptor-based model of type (II. ip exhibiting both pattern 
formation mechanisms, i.e. DDI and hysteresis. We show that the model satisfies the 
theory developed in this paper. Model analysis is supplemented by numerical simulations 
showing de novo formation of far from equilibrium patterns exhibiting jump-discontinuities 
of which the location depends both on the size of diffusion and on the initial data. Section 
[4] presents a detailed proof of the stability theorems stated in Section [2j In Section El we 
prove all statements from Section [3] related to the analysis of DDI initiated formation of 
hysteresis-driven patterns with jump discontinuity. 


2. Main results. Stability conditions for steady states with jump 

DISCONTINUITY 

We state our results on the stability of steady states of (II. ip with finitely many jump 
discontinuities under the following assumptions. 

Assumption 2.1. 

(1) v is a scalar, while u is either a scalar or a vector consisting of two components. 

(2) There exist constants U_, U, V, V, u, u, v, v G M (respectively U_ = (U_i,U_ 2 ), U = 
(Ui, U 2 ), u = («!, u 2 ), u = (ui,u 2 ) G M 2 j such that (i)U_<u<u<U,V_<v< 
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v < V and (ii) 

— oo < u < lim inf min u(t, x) < limsup maxu(i, x) < u < +oo, 

OO x£l £—>-+oo x£l 

— oo < v < lim inf min v(t, x) < lim sup max v(t, x) < v < +oo, 

£—>•+00 x£l £—>■ +OO x£zl 

provided that initial function ( u 0 ,v 0 ) satisfies U_< u 0 < U and V < v 0 < V. The 
inequalities are satisfied componentwise in the case u = 

(3) / and g are twice continuously differentiable on [U_, U] x [V, V] if u is a scalar, 
and on [U_i,Ui\ x [U_ 2 , Uf\ x [V_, V] if u is a vector. 

(4) Initial functions (u 0 ,v 0 ) are in C(I) x (C 2 (/)nC 1 (/)) (respectively C (I) 2 x(C 2 (1)0 
C' 1 (/)), and satisfy Wq(0) = v ' 0 (/) = 0, where I = (0,/). 

Assumption 12. 11 is satisfied if there exists an invariant rectangle, see [241 [25] . 

2.1. The topology. 

In this section, we introduce (eo> A)-stability following the idea of Weinberger in [28] and 
provide conditions for stability of steady states of system (11.11) in this topology. 

For a bounded interval / = (0, /), let BV(I ) denote the space of all functions of bounded 
variation defined on /. We define a neighborhood basis for u 6 BV(I) with values in [u, u], 
where u, u are given by Assumption 12.11 as 

N £A u(u) = {u eBV(7, [u,u]) | 

there exists R C / such that ||w — u\\ 2 Lao ^ < e 2 and meas(J \R) < £ 4 }. 

We are particularly interested in the situation where the initial function u 0 is close to a 
steady state u with finitely many jump-type discontinuities, but u 0 is still continuous on I. 
If we assume v 0 £ C 2 (I) flC(/), then the solution (u(t, x),v(t,x)) of the initial-boundary 
value problem is continuous for all t > 0. If the steady state u is ‘stable’, then u(t,x) 
is expected to converge towards a function with jump discontinuity in space. Therefore, 
the uniform norm is not appropriate to measure the closeness, see Figure 2.1. It turns 
out that N £& u(u) provides us a reasonable topology for this purpose. Following [28], 



Figure 2.1. Illustration of the topology applied to problem (II. ip for scalar 
u. Model (II.ip exhibits steady states with jump discontinuity and global 
existence of classical solutions, u represents a steady state while u(t,x ) 
represents a solution for some t. 


(£ 0 , A)-stability in this topology is defined as 
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Definition 2.2 ((e 0 , H)-stability). A stationary solution (u,v) of system (11.11) is said to 
be (eq, A)-stable for positive constants £o and A if the initial functions (u o,fo) satisfy 

(2.1) ||u 0 - + IN - v\\w (i) < £2 

for some R C I with meas(/ \ R) < e 4 , and for some £ G (0, e 0 ), then 

(2.2) ||«(t, •) - u \\ 1 ^ r ) + IK ■) - v \\ 2 h ^ i ) < Ae2 - 

for all t > 0. Here H l (I) ={w£ L 2 (I) \ u! G L 2 (I)} and 11it11(j-) = || M IU 2 (/) + ||'ii / ||i 2 CD- 

2.2. Stability conditions. 

Definition 12.21 allows us to give conditions for the (e 0 , H)-stability of steady states with 
jump discontinuity. First, we analyze a system of one ordinary differential equation cou¬ 
pled to one reaction-diffusion equation. 


Theorem 2.3 (Stability for a scalar-valued u ). Let I = (0,1) be a bounded interval in M. 
Under Assumption ^. li consider the system of equations 


(2.3) 


du 

for x G I, t > 0, 

dv ^d 2 v . . 

at= D dx 2+sMi 

for x G I, t > 0, 

dv 

Tx (t,x) = 0 

at x = 0, l, t > 0. 


Let (u,v) be a steady state with finitely many discontinuities of u. Denote the Jacobian 
matrix of the kinetic system at the steady state by 


(2.4) 


B(x) 


^^(u( x ),v( x )) 

^(u( x ),v( x )) 


^(u(x),v(x))' 
^{u{x),v(x)) ^ 


and assume that the following inequalities hold for all x G I: 


(2.5) 

(2.6) 
(2.7) 


—(u(x),v(x)) < -Cl < 0 , 
^£(u(x),v(x)) < -Cl < 0 , 


det B(x) > 0, 


where ci is a positive constant independent of x. Then (u,v) is (e q , A)-stable for a pair 
(£ 0) A) with 0 < Eq, A < oo. 

The proof of Theorem 12.31 is similar to that of stability for a particular model with 
cross-diffusion in [28]. We shall present it as a corollary to more general Theorem 12.41 
in Section [4] Theorem 12.31 shows that the steady states of the model are stable, if both 
species are self-inhibitory and cross-inhibition/cross-proliferation is asymmetric. Even if 
each species inhibits all species, a steady state can be stable if the cross-inhibition is 
sufficiently small compared to the self-inhibition. However, in this case, external influx is 
necessary. 












6 


STEFFEN HARTING, ANNA MARCINIAK-CZOCHRA, AND IZUMI TAKAGI 


Next, we consider a model comprising of two ordinary differential equations and one 
reaction-diffusion equation, i.e., u = (ui,u 2 ): 

du i 


( 2 . 8 ) 


dt 

du 2 

dt 

dv 


= fi(u!,U2,v), 

= f 2 (u 1 ,u 2 ,v), 


j d 2 v 

at = D w + 9{UUU2 ’ V ) 

dv dv 

s (*,o) = s ( t ,D = a 


x e /, t > o, 
x E fit > 0, 
x E fit > 0, 


We are interested in the stability question on a given steady state (wi(t), u 2 (x), v(x)) 
of (12. 8 p with finitely many jump discontinuities. It is convenient to put 



a u (x) = 

dfl 

dui 

(ui,u 2 ,v), 

ai 2 {x) 

dfi,~ ~ 

= a (-0 i,U 2 ,V), 

du 2 

a 13 (x) 

dfi,~ ~ ~n 

= -r^{Ul,U 2 ,V), 

(2.9) 

a 2 \{x) = 

df 2 

du\ 

(ui,u 2 ,v), 

022 (x) 

df2 ( ~ ~ ~ x 
= 9m2 (« 

a 23 {x) 

dfo _ 

= -Q X j{ui,u 2 ,v), 


031 (x) = 

dg 

dui 

(ui,u 2 ,v), 

032 (t) 

dg 

= o ("Oi, u 2 , v ), 
du 2 

033 (x) 

dg,~ ~ - x 

= ■^(u 1 ,u 2 ,v), 

and 








(2.10) 




A(x) = 

(Oji (a;) )i<j ) j< 3 . 




Theorem 2.4 (Stability for vector-valued u). Let I C M be a bounded interval. Under 
Assumption \2.1l consider a system of type (12.81) . Let (ui,u 2 ,v) be a steady state with 
finitely many discontinuities of (ui,u 2 ) in x. Denote the Jacobian matrix of the kinetic 
system at the steady state by A{x) as in (12.10]) . Let A^ denote the matrix resulting from 
omitting the i-th row and the j-th column, i.e. Aij = (aki)k^i,i^j- Assume that A{x) and 
Aijfx) satisfy 

3 

(2.11) ti A{x) < 0, tr A(x) det Ajj ( x ) < detH(x) < 0, 

1=1 
3 

(2.12) tr y4 33 (a;) det Ajfix) < detH(x) + tr A(x) det H 33 (:r), 

i=i 

3 

(2.13) 0 < 3 det H 33 (:r) < tr A(x) tr A 33 (x) + det Ajj(x), 

l=i 

(2.14) a 33 (x) < —3k < 0 and detH 33 ((r) > — 3 k tr A 33 (x) > 18k 2 


for all x E I, where k is a positive constant. Then the steady state (ui,u 2 ,v) is (eo,A)- 
stable for some positive constants Eq and A. 


Remark. Condition (12.lip is the Routh-Hurwitz condition for the matrix A(x), i.e., all 
eigenvalues of A{x) have negative real part if and only if (12.lip holds. Note that it implies 

Ej= i det Ajj(x) > 0. 

The above theorem deals with general three-component systems, and it is non-trivial 
to check the set of conditions (12. 1 ID — (12. 141) . practically speaking. However, in the case 
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where one species, say U 2 , reacts very rapidly, the situation becomes more manageable. 
Therefore, we consider the following system with a small positive parameter 6: 


du{ 

~dt 


= fi(u{,u 2 ,v s ) 


for x E I, t > 0, 


(2.15) 


5 


du 2 

~dt 

dv 6 

dt 

dv 6 


= f2(u{,U S 2 ,V S ) 
r) 2 v 


dx 2 


dv 


_(i, o) = _( M ) = ° 


for x E I, t > 0, 
for x E I, t > 0, 
for t > 0. 


Corollary 2.5 (Stability for vector-valued u with a small parameter). Let f = (0,I)cl 
be a bounded interval. Under Assumption 2.1, consider a system of type (12.151b Let 
(ui,U 2 ,v) be a steady state with finitely many discontinuities of ( 111 , 112 ) inx. Let A(x) be 
the matrix defined by (12.1 OH . Let Aij denote the matrix resulting from omitting the i-th 
row and the j-th column, i.e., A^ = (aki)k^i,i^j- Assume that the following inequalities 
are satisfied for all x 6 I: 


(2.16) 

d22(x) < —ki < 0, 

(2.17) 

a3 3 (x) < ~k 2 < 0, 

(2.18) 

det An(x) > 0, 

(2.19) 

det A 33 (x) > -k 3 a 2 2 (x) > 0, 

(2.20) 

det A(x) < 0, 


where k \, fc 2 , fc 3 are positive constants independent of x. Then there exists a positive 
constant 5* such that the steady state (u\,U 2 , v) is (eo. A)-stable for some positive constants 
e 0 and A, provided that 0 < 5 < S*. 


2.3. Quasi-steady state approximation. 

By the quasi-steady state approximation of (I2.15|h we understand system (12.15)) for 5 = 0, 
i.e., the following system of differential-algebraic equations: 


( 2 . 21 ) 


du® 


= fi(ui,u° 2 ,v°), 


dt 

0 = f 2 (u 0 1 ,u 0 2 ,v°), 

^0 _ 

dt ~ dx 2 11 21 h 

dv 0 dv° 

-(t,0) = ^-(t,l) = 0 


x E I, t > 0, 
x G /, t > 0, 

x G /, t > 0, 
for t > 0. 


„ dx v dx 

Assuming that / 2 (m°, u 2 , v°) = 0 has an isolated solution u 2 = u 2 (ufiv°), we can re-write 

dUB as 




du\ 

dt 

dv° 

~dT 

did 

dx 


= *(«?,«;(«?, A A 

dhd 


= D- 


+ g( u i, ul(u 1 , u 2 ), v' 


dx 2 
did 

(t J 0) = —(t,l) = 0 


x E I, t > 0, 
x E fit > 0, 
for t > 0. 


( 2 . 22 ) 
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It follows immediately that any steady state of the reduced system (12.22)1 gives rise to 
a steady state of (12.21)1 . Hence, we want to investigate the following questions: 

(1) Does system (12.1511 exhibit DDI if its quasi-steady state reduction (12.22)1 does? 

(2) Are spatially inhomogeneous steady states of model (12.15P (e 0 , A)-stable if they 
are (e 0 , A)-stable steady states of model (12.2211 ? 

The first question is addressed in Section 4, see Proposition [4761 For df 2 /du 2 < — c < 0, 
it turns out that DDI of the reduced system implies DDI of the unreduced system for 
sufficiently small 5. Here, we address the second question on the ‘transfer’ of stability: 

Theorem 2.6 (Stability transfer). Consider a steady state (ui,u 2 ,v) of the system 

du\ 


(2.23) 


dt 

.dU 5 r, 

I -- 

dt 

dv 5 

dt 


fl(u{,U S 2 ,V 5 ), 

X £ /, t > 0, 

f 2 (u{,u 2 ,v 5 ), 

X £ I, t > 0, 

d 2 v s 

D—+ 9 (uty 2 y), 

x e i, t > o, 


(u{(0,x),u s 2 (0,x),v 5 (0,x)) e C(I ) 2 x (C 2 (/) n <?(/)), 
with homogeneous Neumann boundary conditions for v and its quasi-steady state reduction 
du\ 

=D gn + gK , u . Ki AA x g i,t > o. 


=fi(u°,u* 2 (u° 1 ,v°),v 0 ), 
d 2 v° 


x e /, t > o, 


dt dx 2 

Denote the Jacobian matrix of the kinetic system of (12.2311 with 5 = 1 at (ui,u 2 ,i ’) by 
A = (a,ij) and denote the Jacobian matrix of the kinetic system of ()2.24j) at (ui,v) by 
B = (bij). Suppose that 

(1) system (12.24)1 satisfies the conditions of Theorem ] 2. A 

(2) system (12.23)1 satisfies Assumption ] 2. 1[ 

(3) <222, «33 < —c < 0, 

where c is a positive constant independent of x. Then, there exists a 5* > 0 such that 
(ui,v*(ui,v),v) is an (e 0 , A)-stable steady state of system (I2.23P for all S < 5*. 

3. Application to a system exhibiting DDI and hysteresis 

In this section, we present a model exhibiting the coexistence of DDI and hysteresis. 
First, we propose a model consisting of two ordinary differential equations coupled to one 
reaction-diffusion equation. We reduce this model using quasi-steady state approximation. 
Then, we rescale the reduced model and show that it satisfies the conditions of Theorem 
12.31 Finally, based on Theorem 12.61 we show that the stability of steady states of the 
reduced model implies stability of corresponding patterns in the original model. 

Our example takes the form 


= — v\u\ — j3u\v 6 + 9 1 


u 


S\2 


(3.1) 


du[ 

~di 

.dU S 0 

i— 4- 
dt 

dlJ ^ LJ U 

—F = 7 ytt - YsT 5 - fdu[v 5 + 6*2 


1 + k(u\) 2 
= —u 2 u 5 2 + fdu\v 5 — au 2 , 
d 2 v s 


+ au 2 , 


x e I, t > o, 
x e i,t > o, 


<? 


dt 


dx 2 


1 + 2c(nf) 2 


■f" OLVb^-, 


x G /, t > 0, 
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supplemented with homogeneous Neumann boundary conditions and initial functions 
(uf (0), 7 / 2 ( 0 ), v s (0)) G C(I ) 2 x (C 2 (J) flC(7)). Here, a, (3, 7 , 9\, 9 2 , zq, u 2 , v 3 are positive 
constants. Without loss of generality, we consider / = (0,1) in order to simplify analysis. 

The model is a modification of the receptor-based model proposed in [121 US]- Nu¬ 
merical investigations for k = 0 show spike solutions blowing up in infinite time, see [7]. 
Introducing k > 0 reflects saturation of de novo production. The modification assures 
uniform boundedness of the model solutions. 

Variables U\ and u 2 describe cell surface receptors at free and occupied states, re¬ 
spectively. The state of receptors is being changed through binding to diffusive ligands, 
denoted by v, and a reverse process of dissociation of bound receptors. Free receptors and 
ligands are produced de novo, what is described by a function with saturation effect for 
k > 0 . 

For h = 0, (13. ip reduces to a system of differential-algebraic equations. In this case, the 
second equation of (13.ip can be solved uniquely for leading to a system of one ordinary 
differential equation coupled to one reaction-diffusion equation. This reduction leads to 
simplification of analytic investigation of the model. We show that the limit 6 —» 0 is 
regular in the sense of stability-preservation of steady states. 

The quasi-steady state approximation, i.e. 5 = 0, reads 


(3.2) 


dui 

~di 

dv° 


"2 


0\2 


v 2 + a 


= -v 1 u\ - 13 

d 2 v° 0 
= - OT - f3 


u y 


+ 9 1 


K) 


1 + k(m°) 2 ’ 


v2 


0„,0 


dt ' dx 2 
After rescaling u\,v°,t, we obtain 

du 
dt 


v 2 + a 


(3.3) 

(3.4) 


= —u 


uv + m 1 


u\v 


u 


+ 9-2 


K) 5 


1 + k(UiY 


1 + ku 2 ’ 


dv 

dt 


= D 


d 2 v 

dx 2 


/i 3 u — uv + m 2 


u 


1 + ku 2 ’ 


x E I, t > 0, 


x G I, t > 0, 


x G I, t > 0, 


x G I, t > 0, 


(3.5) 

(3.6) 


dv dv 

0(0, x), u(0, x)) = (u 0 (x),v 0 (x)) G C(I) x (C 2 (I) n C(I)), 


with positive initial data. 

To simplify notation, we define 

u 2 

fr(u,v) ■■= -u-uv + m l i + ku2 , 


g r (u, v) := ~n 3 v uv + m 2 ^ ^ 2 . 

We show that the model has a spatially homogeneous steady state that is stable to spatially 
homogeneous perturbation and is destabilized due to autocatalysis in dynamics of non- 
diffusive component. Due to the nonlinear nature of the reaction terms, there exists a 
‘far from equilibrium’ regime, where the species become self-inhibitory, stabilizing steady 
states with jump discontinuity. 

Our first step is to show that models (I3.3[) - (I3.5D and (13.1[) satisfy Assumption 12.11 
i.e., the solutions are uniformly bounded for nonnegative initial functions. Numerical 
simulations of the model solutions are presented in Figure EZ3 We formulate the results 
in the remainder of this section and defer their proofs to Section 0 
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3.1. Global existence and uniform boundedness. 

For model (13.3|) - (13.61) . global existence of a solution in G 1 ((0, oo); C(I ) x ( C 2 (I ) fl C(I ))) 
is assured within the framework of invariant rectangles, see, e.g., ra. The solution is 
uniformly bounded, i.e., 

Lemma 3.1 (Uniform boundednes of solutions of system (I3.3l) - (l3.5p h (i) Assume that 
mi > 2 \[k. Then, for positive initial functions, the set 

(°4(».i+vUF«))x(o^) 

is an attracting invariant rectangle of model ra-flUD. (h) Assume that m\ < e l\fk. For 
initial functions (uo(x), vo(x)) satisfying 0 < uo(x), 0 < i’o(x) < m 2 /(kp 3 ), it holds that 
u(t,x ) is monotone decreasing in t E (0, oo) and 

(3.7) (u(t,x),v(t,x)) —> (0,0) uniformly on I, as t —>■ +oo. 

Also for model (13.ip . global existence of a solution in C' 1 ((0, oo); C(I) 2 x (C 2 (I)nC(I))) 
follows from the theory of invariant rectangles. 

Lemma 3.2 (Uniform boundednes of solutions of system (13. ip ). For positive initial func¬ 
tions, the set 



x 0, 


e i 


kS min(zq, u 2 ) 


x 




aOi \\ 
5a min(zq, u 2 )) ) 


is an attracting invariant rectangle of model (ED. Here, in the case Q\ < 2 vx^/k, the 
expression 6 l 1 /( 2v{) + ^ (6 l 1 /(2z/ 1 )) 2 — k is interpreted as an arbitrary positive number. 
For 9 1 < it holds that 

{u{(t, x), u s 2 (t, x), v 5 (t, x)) —^> (0, 0, 0) uniformly on I, as t —> +oo. 


3.2. Existence of steady states. 

By definition of quasi-stationary reduction and linear rescaling, there exists a one-to- 
one mapping between the sets of steady states of model (I3.3l) - (l3.5p and of model (13.11) . 
Therefore, we may focus on existence of steady states of model (I3.3p ~ (13.51) . A function 
( u,v) is a steady state of model (13.3p ~ (13.51) only if 

du 

(3.8) _ = /r(M , v)= 0, 

is satisfied. Note that u{v) := {u G M>o | f r (u,v) = 0} may consist of more than one 
element. We characterize different possible branches of solutions for given w. 


Lemma 3.3. For given v G R> 0 ; 

(3.9) 


f r {u,v) = 0 
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has three solutions 

(3.10) u 0 (v) = 0, 

For v > mi/(2\/k) — 1, U-(v) and u+(v) are not real-valued. We therefore define 
(3 - 13) 

which is the largest v G M such that u±(v) is real-valued. Note that v r > 0 is equivalent 
to mi > < l\fk. Therefore, v r < 0 implies (u(t,x),v(t,x)) —> (0,0) uniformly as t —> oo, 
due to Lemma P 

We show in Lemma 13.71 that the steady states have different stability properties de¬ 
pending on the branch. 

Definition 3.4. We say that a steady state (u,v) is of type U- (respectively u + , uq) at 
x G / if u(x) = U-(v(x)) (respectively u(x) = u + (v(x)), u(x) = uq(v{;x)) = 0). 

The existence of steady states can be summarized in the following lemma, illustrated 
in Figure EID 



0 

0 


Figure 3.1. Plot of the nullclines of f r (u,v) = — (l+v)u+mi(u 2 /(l+ku 2 )) 
for m ^ 0 and g r (u, v ) = —(/i 3 + u)v + m 2 {u 2 / (1 + ku 2 )). 

Lemma 3.5 (Existence of spatially homogeneous steady states). 

(1) For arbitrary positive parameters, the trivial solution (u,v) = (0,0) is a spatially ho¬ 
mogeneous steady state of model (13.31) - (13.51) . 

(2) Let mi < m 2 and 



Then, there exists a positive k\ such that for all k < k*, system (I3.3[) - (I3.5[) has 
exactly two positive homogeneous steady states. Moreover, if k < min(/;;(, ((m 2 — 
mi)/(mi/i 3 )) 2 ), both nontrivial positive homogeneous steady states are of type ?j_. 
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We show in Lemma 13.71 that the conditions of Theorem 12.31 are never satisfied if u is of 
type U- at some x G /, but always satisfied if u is of type uq or u+ for all x £ I. Hence, 
the latter type of solutions is constructed. 

Lemma 3.6 (Existence of steady states with jump discontinuities). Assume < l\[k < mi < 
m 2 . Then, there exist infinitely many steady states (u,v) G (L°°(J) n BV(I )) x ( H 2 (I ) D 
C 1 (/)) which are of type u + or uq for all x G /. 


3.3. Stability of steady states for the quasi-steady state approximation. 

Following the same principle as in Section [2j we prove stability of steady states of the 
quasi-steady state approximation fl3.3H - fl3.5H . In the next subsection, we conclude stability 
of steady states of model (13.11) . In order to apply Theorem 12.31 we investigate how the 
entries of Jacobian matrix at a steady state depend on the branch: 


Lemma 3.7 (Characterization of Jacobian matrix for model fj3.3j) - fl3.5fl ). Assume 2 \[k < 
mi < m 2 . For a pair of nonnegative functions (u(x),v(x)), denote the Jacobian matrix 
of the kinetic system of (I3.3p - (l3.4p by 


B(x) = B(u(x),v(x)) = & 


J9r , 


39r , 


x, 


< 3 - 14 > . _ 

du V X 1 dv 

Then, it holds for 0 < v(x) < mi/ (2 y/k) — 1 v r that 


ZF( x ) (x) 


bu(x) hfix)' 
b 2 i{x) b 22 (x) I 


(3.15) 


bzi(x) 


> 

0, 

lfu{ 

[x) 

= U-{v( 

») 

1 

r < 

o, 

ifu[ 

» 

= u_(v 1 

^)) 

< 

0, 

ifu{ 

[x) 

= u + (v{ 


bi2(x) i 

< 

o, 

ifu{ 

>) 

= u + (v 1 

{*)) 

< 

0, 

if u{ 

>) 

= Uq 


1 

1 = 

0, 

if u[ 

[x) 

= Uq 



0, 

o, 



/ . 

1 

1 

r < 

o, 

if u{ 

\x) 

= u_(v 1 

(x)) 

> 

< 

if u{ 
if u[ 

a) 

\x) 

= u_(v 1 

= U 0 

kni?) < 

< 

0, 

if u{ 

» 

( \ 

= u + (v 1 

{%)) 






l < 

o, 

Vu 1 

a) 

= Uq 



Moreover, there exists 0 < v* < v r such that, for u(x) = u + (v(x)), it holds that 
(3 16) h (r) I >0 E(v*,v r ), 

(3 ' 16) b 2 i{x)^ Q zfv ^ = v *' 

For an isolated solution u(v) of f r (u,v) = 0 (i.e., df r /du ■ df r /dv 0 at (u(v),v)), it 
holds that 


(3.17) 


det(B(u(v),v)) = —b 12 ( fu(v) 


-1 


d 

dv 


g r (u(v),v), 


and, in particular, for u(x) = u + (v(x)), 

(3.18) det(H(M + (f), v)) > 0. 


Applying Theorem 12.31 yields stability of the steady states constructed in Lemma 13.61 


Corollary 3.8 (Stability of solutions with jump discontinuity of model (I3.3lf - (I3.5IL . 
Assume D > 0 and 2 \fk < mi < m 2 . Then, there exist infinitely many (e 0 , A)-stable 
steady states with jump discontinuity of model fl3.3p - fl3.5ll . which are for every x G / 
either of type u 0 or u + . 


To obtain DDI, it is left to prove that at least one of the spatially homogeneous steady 
states of type w_ is stable under spatially constant perturbations, i.e. a stable steady state 
of the kinetic system. Instability of all spatially homogeneous steady states of type w_ 
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follows from |T6j due to autocatalysis of the non-diffusive species. Existence of exactly 
two homogeneous steady states of type w_ results from Lemma 13.51 

Lemma 3.9 (Diffusion-Driven Instability in model (I3.3p - fl3.5p ). Let mi < min(m 2 , y/mjj)- 
If there exist two homogeneous steady states (u-(v-),v-), (u-(v + ),v + ) withv- < v + , then 
there exists k 2 > 0 such that for all k < k 2 

(1) (u-(v-),v-) is unstable to spatially homogeneous perturbations, 

(2) (u-(v + ), v + ) is stable to spatially homogeneous perturbations. 

Furthermore, all spatially homogeneous steady states of type U- are unstable for k > 0 
and D > 0. 

We summarize our results in the following 

Theorem 3.10 (Coexistence of DDI and Hysteresis for model (I3.3p - fj3.5p ). Let k\ and k 2 
be defined as in Lemmas \3.5\ and \3.9\ Under conditions 


(3.19) 


m 1 ,m 2 ,k,fi 3 > 0, 

mi < min(m2, 



there exists a k* in the interval 


(3.20) 


0 < k* < min I k*, k 2 , 


m 2 — mi 
mi/i 3 


m t 


such that, for all 0 < k < k*, the following (l)-(5) hold true: 

(1) system (j3.3p - (13.5p has exactly two strictly positive spatially homogeneous steady 
states (w_(u_),u_) and (u_(v + ),v + ) with u_ < v + . 

(2) (w_(u_), u_) is an unstable equilibrium of the kinetic system for fl3.3p - fl3.5p . hence 
is an unstable steady state of system (13.31) - (13.51) . 

(3) (w_(i> + ),u+) is a stable steady state of the kinetic system of fl3.3p - fl3.5p and an 
unstable steady state of (13.3D - 03.51) . 

(4) (0, 0) is a stable steady state of the kinetic system for (I3.3p - fl3.5p and of the original 
system (lT3D-(im 

(5) system fl3.3p ~ fl3.5p has infinitely many (e 0 , A)-stable steady states with jump discon¬ 
tinuity, which are at x G / of type u + or uq . 

A numerically obtained solution, using the finite element library deal.ii |0], is presented 
in Figure [3~2l We observe destabilization of the spatially homogeneous steady state and 
convergence towards a pattern with jump discontinuity. Compartment u is discretized 
with cell-wise constant finite elements and v is discretized with cell-wise linear, globally 
continuous finite elements. The formulation is the weak formulation. The time-stepping 
scheme is Crank-Nicholson. 
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Figure 3.2. Numerically obtained solution to model (I3.3H - (13.5|) for 
parameters m\ = 1.44, m 2 = 2 , /i 3 ps 4.1, k = 0.01,i4 = 1. We 
observe convergence towards a steady state with jump discontinuity. 
Left: Non-diffusive component u. Right: Diffusive component v. 


3.4. Stability of steady states for the unreduced model. 

In this section we show that model (13.111 satisfies the conditions of the ‘stability transfer’ 
Theorem 12.61 


Lemma 3.11. Consider model (13. ip . Define U 2 (ui,v) := fiuiv / (u 2 +a). Under conditions 

(3.21) 0 < z/i, 1/3,01, 02,^2, /?,7> 

(3.22) 9i < min ( 0 2 , x /9 2 fi- 


V2 


(3.23) 


z / 3 > 2vif3 


a + v 2 


v 2 1 I 26o — 9\ 


a + z / 2 9i 


0i 


+ 2W -F 


there exists a constant k* in the interval 


(3.24) 


0 < n* < min ( ft 


, 6*2 - 9 1 z / 2 

ui a + z / 2 


9i_ 

2ui 


and a positive constant 5* such that for all 0 < 5 < 5* and for all 0 < n < k* the following 
statements hold true: 

(1) system (I3.1|) has exactly two strictly positive spatially homogeneous steady states, 
(u-(v-),u 2 (u-(v-),v-),v-) and (u-(v + ),U 2 (u-(v + ),v + ),v + ) with V- <v + . 

(2) (m_(u_), U 2 (u-(vJ), u_), u_) is an unstable steady state of the kinetic system for 
(13.ip and of the original system (13.11) . 

(3) (u_(u + ), U 2 (u_(v + ),v + ), v + ) is a stable steady state of the kinetic system of (13.11) 
and an unstable steady state of (13.ip . 

(4) (0, 0, 0) is a stable steady state of the kinetic system for (13.11) and of the original 
system (|23J . 

(5) system (13.ID has infinitely many (£ 0 , A)-stable steady states with jump discontinuity, 
with component Ui being at x E I of type u + or u 0 . 
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3.5. Numerical results: Effects of initial conditions and diffusion coefficient. 

In the previous subsection, we showed that a spatially homogeneous steady state is desta¬ 
bilized and that infinitely many, discontinuous steady states are stable, both due to in¬ 
troduction of diffusion. In this subsection, we show some numerical results showing that 
the arising pattern depends on both, the initial conditions and the diffusion coefficient. 
In figure PTTTH numerical approximation of solutions to problem (I3.3l) - fl3.5l) are shown. We 
observe that for the same initial conditions and parameters of the kinetic system, the 
number of jump-type discontinuities varies according to the size of diffusion coefficient D. 
A trend towards a higher number for smaller diffusion-coefficient can be observed and the 
shape of the arising pattern is different from the shape of the initial conditions. For large 
diffusion coefficient, a ‘plateau’ arises close to the spatial position on which the initial 
condition of u assumes its maximal value. Here, initial conditions seem to determine the 
position of ‘plateaus’ as can be observed in figure 13.41 Note that a local maximum is 
positioned at x = 0 in figure 13.41 



Figure 3.3. Numerically obtained solution component u of model 
(I3.3p - 03.5I) for parameters rn i = 1.44, m 2 = 2, /i 3 ~ 4.1, k = 0.01 with 
varying diffusion coefficient: upper left: D = 5, upper right: D = 1, 
lower left: D = 0.5, lower right: D = 0.1. We observe emergence of 
more jump-type discontinuities for smaller diffusion coefficient. Ini¬ 
tial conditions are u 0 (x ) = 1.725 — 0.1 cos(27ra: 2 ), v 0 (x) = 2.48615. 
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20 


0 


0 



10 o 


1 


Figure 3.4. Numerically obtained solution component u of model 
fl3.3p - fl3.5p for parameters rni = 1.44, m 2 = 2, p 3 cs 4.1, k = 0.01, D = 
5. We observe emergence of jump-type discontinuities around lo¬ 
cal maxima of the initial conditions. Initial conditions are Uo(t) = 
1.725 - 0.l£ 4 cos(87T£ 2 ),u 0 (>) = 2.48615. 


4. Proof of the stability theorems 


4.1. Spectral structure of the linearized operator. 

To handle Theorem 12.31 we introduce the following notation for simplicity: 



Lemma 4.1. Under the assumptions of Theorem \2.3\. consider the operator 



(4.1) 


on (L 2 (/)) 2 with domain L 2 (I) x ( H 2 (I ) D {f> \ if' = 0 at x G dl}). Then —Ci is a 
sectorial operator on (L 2 (/)) 2 , i.e., there exist positive constants M 1; n 0 and oj e (0,7t/2) 
such that 


(4.2) || (A + Cf) 1 |l ij( ( L 2 (J) ) 2 ) < i A _^ o j +ft . o f° r all A G {A G C | |arg (A — k 0 )| > ca}, 


where B((L 2 (I)) 2 ) denotes the Banach space of all bounded linear operators on (L 2 (/)) 2 
equipped with the operator norm. Hence C\ generates the analytic semigroup e tCl on 
(L 2 (/)) 2 and satisfies the estimate 


(4.3) 



for all t > 0. 


Here M and are positive constants independent of 

Lemma 4.2. Under the assumptions of Theorem\2.f\ consider the operator 
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defined on (L 2 (/)) 3 with domain (L 2 (/)) 2 x ( H 2 (I ) D {if \ 0' = 0 at x G dl}). Then 
—C -2 is a sectorial operator on ( L 2 (I )) 3 , i.e., there exist positive constants M 2 , k 0 and 
u> G ( 0 , 7 t / 2 ) such that 

i 11 3 Mo 

(4.5) ||(A + C 2 ) L( (L 2 (/)) ) < |^_ Ko | + /to f° r al1 A e ( A e C I |arg(A- k 0 )| > w}, 

where £>((L 2 (/)) 3 ) denotes the Banach space of all bounded linear operators on (L 2 (/)) 3 
equipped with the operator norm. Hence C 2 generates the analytic semigroup e tc2 on 
(L 2 (/)) 3 and satisfies the estimate 


(4.6) 



< Me Klt (||0i||l2 (J ) + ||0 2 || L 2 (J) + U\\lHd) 

(L 2 (I )) 3 


Here M and Ki are positive constants independent of (0i, 0 2 ,0). 


for all t > 0. 


Since Lemma 14.11 is obtained as a corollary to Lemma 14.21 we prove only Lemma 14.21 
(Indeed, choosing a > 0 sufficiently large, we see that the matrix 


f-a 0 0 A 

0 b n (x) b 12 (x) 

\ 0 & 2 i(ic) b 22 (x)J 


satisfies conditions (12. 1 if) (12. 14j) of Theorem 12.41 1 

Proof of Lemma Instead of —C 2 we study C 2l the estimate (14.5)) is obtained by 
replacing A with —A in il 1. 1 III below. Hence, we consider the nonhomogeneous equation 


(4.7) C 2 (t 2 j = A ^j 

for r±, r 2 , s G L 2 (/) under homogeneous Neumann boundary conditions on 0. 

Step 1. By assumptions (12.11 jl — (12. 13[) we see that trH 33 (a;) < 0 < detH 33 (x), and 
hence both eigenvalues pfix) and p 2 (x) of H 33 (x) lie in the left-half plane. In what fol¬ 
lows we number p 1; p 2 so that Re /ifix) < Re p 2 (x). More precisely, there exists a positive 
number p* such that the rectangle 1Z = {A G C | — p* < Re A < — 3 k, |ImA| < /i*} con¬ 
tains pi(x) and p 2 (x) for all x G I. Indeed, pi fix) are bounded because aij(x ) are bounded. 
It is elementary to check that R ept,j{x) < maxjtr H 33 (x)/2, det H 33 (x)/tr A 33 (x)}. By as¬ 
sumption (12.141) . the right-hand side of this inequality does not exceed —3 k. Therefore, 
in particular, 

<t(H 33 (o;)) c77c{AGC|ReA< -3 k}. 

Hence, whenever Re A > —3 k, we can solve 


(4.8) 


for (0i ,(j> 2 ): 


(4.9) 


01 = 


02 — 


(H 33 (a;) - A) (^j 




det 


det 


fr\ a 12 (x) A 
\r 2 a 22 (x)-\J 

fan(x) — A rfi\ 
V 021 00) r 2 J 


— det 

— det 


fa i 3 (x) a 12 (x) A 

0023 00) a 22 (x) — Ay 

fau(x)-\ oi 3 0 c)A 
V «21 (x) 023 (x)J 



det(H 33 (a:) - A), 



/ 


det(H 33 (a;) - A). 
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Inserting this result into the third equation yields a scalar inhomogeneous Sturm-Liouville 
problem subject to homogeneous Neumann boundary conditions: 


(4.10) 
where 

(4.11) 


D 


d 2 ip 
dx 2 


q(x, X)ip = pi(x, A)r x + p 2 (x, X)r 2 + s, 


det Ai 3 (x) + Xa 3 i(x) det^4 23 (^) - Aa 32 (x) 

pi{x,A) — - , .—7T-rr—, P 2 \x,A) — 


q(x, x) = 


det(^(x) - A) 
det(44(x) — A) 


det(^(x) - A) 


det(.A 33 (x) - A)' 

Step 2. Let us define a sesquilinear form B(ip i,ip 2 ) on H l (I ) by 


(4.12) 




dip 1 dip 2 
dx dx 


dx + / g(x, A)^ ip 2 dx. 


Sublemma 4.3. There exist positive constants and T\ such that if A G C satisfies 
Re A > —2 k 0 or |Im A| > t/ien i/ie sesquilinear form B is bounded and coercive. 
Namely, there exist positive constants M 0 , y 0 such that 


(4.13) \B(i> 1 ,ip 2 )\ < M 0 ( 1 + |A|) \\ipi\\m(i)\\ifa IIhrj) / or all ipi, ip 2 e H X (I), 

(4.14) ReB(ip,ip) > 7 o||V , ||^i(/ ) for all ip e H 1 (I). 

Here, 70 depends on X. 


Proof of Sublemma \f.3\ First we prove the boundedness of B. 
introduce the following notation: 


It is convenient to 


(4.15) 


" 3 

D(x) = det 44 (x), T(x)=trA(x), S(x) = det 44 j 3 -(x), 

3 =1 

L >3 (x) = det 44 33 (x), T 3 (x) = tr 44 33 (x), 

k P( A) = det(A(x) - A), Q(A) = det( 44 33 (x) - A). 


We note that P{ A) = Z)(x) — 5(x)A + T(x)A 2 — A 3 and Q(A) = D 3 (x) — T 3 (x)A + A 2 , so 
that 


(4.16) |g(x, A)| = |P(A)/Q(A)| = |A| + 0(1) as |A| ->■ +00 

uniformly in x € /. Moreover, zeros of Q(A) are both confined in the rectangle TZ defined 
in Step 1 for all x G I. Therefore, there exist positive constants do an d itiq such that 

|g(x, A)| < m 0 (|A| + 1) 

for all A satisfying Re A + 2k, > 0 or |ImA| > To, which implies the desired inequality: 
I^C0i>^ 2)1 < ^ > IIV ,/ iIIl 2 (7)||'02IU 2 (i) + m o(|A| + i)||V’i||l 2 (j)||V’2||l 2 (7)- 
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Second, we prove the coerciveness of B. Let us start with the case where |A| is sufficiently 
large. From 

(4.17) P(A)Q(A) + P(A)Q(A) = 2D(x)D 3 (x ) - (D 3 (x)S(x) + D(x)T 3 (x))(X + A) 

+ {D 3 {x)T(x) + D(x))( A 2 + A 2 ) + 2S'(x)r 3 (a;)|A| 2 
- D 3 (x)( A 3 + A 3 ) - (T(x)T 3 (x) + 5'(a;))|A| 2 (A + A) 

+ E 3 (;r)|A| 2 (A 2 + A ) + 2T(x)\A\ i — |A| 4 (A + A), 

it follows, due to A 2 + A 2 = (A + A) 2 — 2|A| 2 and T(x) — T 3 {x) = a 33 (x), that 

(4.18) -ReP(A)Q(A) = |A| 4 (Re A - a 33 (a;) - 2T 3 (a;)(Re A) 2 /|A| 2 + OdAI" 1 )) 
as |A| —> +oo, uniformly in x G I. On the other hand, we have 

(4.19) |Q(A)| 2 = D 3 (x) 2 - 2 D 3 (x)T 3 (x) + ( T 3 {x ) 2 - 2D 3 (x))\ A| 2 

+ 4D 3 (x)(Re A) 2 - 2T 3 (a:)|A| 2 Re A + |A| 4 . 

Therefore, we see that 

(4.20) Re q(x, A) = Re A — a 33 (x) — 2T 3 (x)(Re A) 2 /1A| 2 + 0(|A| _1 ) as |A| ->• +oo, 

uniformly in x e I. Since Re A — a 33 (x) > Re A + 3 k and T 3 (x) < 0 by assumption (12.141) . 
we conclude that there exists a positive constant /\ > Jq such that 

(4.21) R + (Re A + 3 k)||^|| 2 2 (/) > min{D, k/ 2}\\^\\ 2 h1{i) 

whenever Re A > —2k and |ImA| > J\. 

Next we turn to the case where |ImA| < J\. By direct computations we have 

Re det(R(x) — A) = det(y4(x) — Re A) + (ImA) 2 (3ReA — T(x)), 

Im det(A(a;) — A) = (ImA)((ImA) 2 + (2 T(x) — 3ReA)ReA — S(x)), 

Re det(R 33 (a:) — A) = det(A 33 (x) — Re A) — (ImA) 2 , 

Im det(y4 33 (x) — A) = (ImA)(2ReA — T 3 {x)). 

Using these expressions, we obtain 

( 4 . 22 ) for A = € + i„, 

where we have dehned as follows: 

Qsiv) = !, 

Qa(v) = ~ T ( X ) ~ T s( x ), 
q 3 (v) = D s( x ) + S(x) + T(x)T 3 (x) + 2 ?y 2 , 
q2(v) = ~ { D (x) + D 3 {x)T{x) + S{x)T 3 (x) + 2 T(x)?f) , 
qM = D 3 (x)S(x ) + D(x)T 3 (x) + (T(x)T 3 (x) + S(x) - 3 D 3 (x))rj 2 + rf, 
q 0 (v) = ~D(x)D 3 (x) + (D(x) + D 3 (x)T{x) - S(x)T 3 (x))rj 2 - (T(x) - T^x))^. 
Thanks to assumptions (12.lip — H2.14j) of Theorem 12.41 we see that 

(4.23) q 3 {rj) = 1, q 4 (v) > 0, q 3 (v) > 0, q 2 (v) > 0, q^rj) > 0, q 0 (rj) > 0 
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for all x E I and Re A + 2k > 0. Observe that 


(4.24) ^ + q4 ( v )e = ^-(T + T 3 )y, 

(4.25) q 3 ( V )e + q 2 ( V )e > £ 2 ((r 3 + TT 3 + S + 2 rf)^ - (D + D 3 T + ST 3 )); 

(4.26) qi (r,)f + qo(ri) > {D 3 S + DT 3 + (TT 3 + S- 3D 3 )rj 2 + r? 4 } £ - DD 3 . 


Let n 0 be a positive constant defined by 

O f trrr \ , w D (x) + D 3 (x)T(x) + S(x)T 3 (x) 

I xei Ji D 3 (x) + T(x)T 3 (x) + S{x) + 2 T{ 


2 ’ 


sup 

16/ 


D(x)D 3 (x) 


D 3 (x)S{x) + D(x)T 3 (x) + (T(x)T 3 (x) + S(x) - 3 D 3 ( x ))r 2 + r 4 


We may assume that kq < k. If £ > —3kq and |ImA| < A, then 


£ — {T + T 3 ) > £ + 3k 0 ; 

(D 3 + TT 3 + S + 2if)i — (D + D 3 T + ST 3 ) 

> (A. + rr 3 + 5 + 2// 2 )e + 3 a 0 {D 3 + TT 3 + S + 2T 2 ) > {D 3 + TT 3 + S)(£ + 3a 0 ); 
{(D 3 S + DT 3 ) + (TT 3 + S' - 3 D 3 ) V 2 + r? 4 } e - DA 

> { (D 3 S + DT 3 ) + (TT 3 + 5 - 3D 3 )g 2 + if} £ 

+ 3a 0 { (D 3 S + dt 3 ) + (TT 3 + r- 3 d 3 )a 2 + rf} > (d 3 s + da)(£ + 3« 0 ). 

Therefore, 

£ 5 + qfv )£ 4 + <?s(h )£ 3 + Q2(v )£ 2 + <?i(h)£ + Qo(v) 

(4 27) 

> (£ + 3« 0 )(e 4 + (A(x) + T(x)T 3 (x) + S(x))£ 2 + D 3 (x)S(x) + D(x)T 3 (x)). 


Here we note the following elementary lemma, the proof of which is omitted. 


Lemma 4.4. Let /(£) = £ 4 + a 3 £ 3 + a 2 tf + ai£ + ao and g(£) = 5(£ 4 + f3 A 2 + fif) be two 
polynomials, where /3 q and fii are positive. Then there exists a positive constant b 0 such 
that /(£) < g(£) for all f ifb>b 0 . 

By applying this lemma to |Q(A)| 2 and £ 4 + (D 3 + TT 3 + S)f 2 + D 3 S + DT 3 , we see 
that 


|Q(A)| 2 < C^ 4 + (D 3 (x) + T(x)T 3 (x) + r(x))^ 2 + D 3 (x)S(x) + D(x)T 3 (x)) 
for some positive constant C\. Consequently we obtain 

(4.28) Re q(x, A) > ^ + ^° 

t'i 

provided that Re A > —3 k 0 and |ImA| < A. Hence, 

(4.29) Re^,^) >^11^11^(7)+ ^r 4 (ReA+ 3/co)||^||i 2( r) >min{D,«o/Ci}||^||^ (J ) 

whenever Re A > — 2kq and |ImA| < A. Combining (14.29)1 with (I4.21j) completes the 
proof of the coerciveness estimate (14.141) . hence Sublemma 14.31 □ 

Step 3. By the Lax-Milgram theorem, for any h G L 2 (I) there exists a unique A. £ 
H l (I ) such that 

(4.30) 


= (AVOiR/) for a11 ^ e 
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as long as Re A + 2k 0 > 0 or |ImA| > A- We would like to prove that there exists a 
positive constant M\ such that 

whenever Re A > — k 0 . 

By (FOH with -0 = i/) h , (POD and (0Z£1 we obtain 

(4.31) D\\^' h \\ 2 L 2 {I) +7o(ReA + 2K 0 )||A||i2 (/) < Re B(ip h , ip h ) < ||h|| L 2 (7) ||0 fe || L 2 (/) . 

As in the computation of (14.17[) . we see that 

(4.32) 

ImP(A)Q(A) = -(ImA) {(D 3 (A)S(A - D(x)T 3 (x)) + (D(x) - D 3 (x)T(x)( A + A) 

-D 3 (x)(X 3 - A 3 ) - (T(x)T 3 (x) - S0r))|A| 2 - T 3 (x)|A| 2 (A + A) + |A| 4 } , 

which yields 

(4.33) Im q(x, A) = Im A(1 + 0(|A| -1 )) as |A| —» +oo. 

Hence there exists a positive constant A(> A) such that if |ImA| > A then 


|ImH(A,0/i)| = 


Im q(x, A) 10ft| 2 dx 


> 


| Im A | 


Uh\\l\iy 


Combining this estimate with (14.301) . we conclude that 

( 4 -34) |ImA| U h f L2{I) < 2||h|| L 2 (7) ||^|U2 ( r ) . 

From fl4.3ip and (I4.34p we immediately obtain, as long as |ImA| > r 2 , 

(4.35) ( 7 o 2 (Re A + 2 % ) 2 + (Im A) 2 )||*||i 1(/) < 5||ft||l =(J) . 

If Re A + Ko > 0, then (Re A + 2 k 0 ) 2 > (Re A + kq) 2 + Kq. Hence, 04.35p implies 

mi 


(4.36) 


||Ai||l2(/) < 




|A + Kq| + Kq 

for Re A + k 0 > 0 and |Im A| > A- Then from 04.31 p and (02 

(4.37) mu//ii2 - ,? 'i 


it follows that 


D 1I0 /iIIl2 (/) < ix , 

v ' |A + K 0 \ + K 0 


2 

L 2 (iy 


Finally we treat the case |Im A| < A and Re A + kq > 0. Then we have Re A + 2kq > 
7 i(|A + Ko| + ^o) with 2^\ = k\/[T 2 + Kq). Hence, 04.31 p yields the desired estimate 

(4.38) T>||V^|| 2 2 (7) + 7o7i(|A + k 0 | + no)\\^h\\ 2 L 2 {I) < ||h|| L 2 (J) ||'0 ft || £ ,2 (J) . 

Step 4- Now we are ready to finish the proof of Lemma 14.21 Clearly from 04.lip it 
follows that 


(4.39) sup |pi(a:, A)| < C 2 and sup | p 2 (x, A)| <C 2 

ccE/ x£l 

for some positive constant C 2 independent of A ^ 77. Therefore, the right-hand side of 
(| 4.10P is bounded by A 3 (||ri|| L2(7) + ||r- 2 || L 2 (J) ) + ||s||l 2 (/)- 

Hence, the unique solution -0 of Q4.10p under homogeneous Neumann boundary condi¬ 
tions satisfies the estimate 

(4.40) || 0 ||l 2 (/) < lx , M \, — {CsdlnW^ii) + INU 2 (.o) + IM|l 2 ( j)}. 

| A K,q | Kq 
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From (14.911 we see that 

||0l|U 2 (J) < Cn(A)||r 1 || L 2 ( 7 ) + C i2 (A) ||r 2 ||z, 2 (T) + Ci 3 (A) 1111 Z .2 (j) , 
||02||l 2(/) < c 2 i(A)||ri|| L 2 ( 7 ) + c 2 2 (A)||r 2 ||z, 2 (t) + c 2 3(A)1111z, 2 (j), 


where 


cn(A) = 
Cl3(A) = 
C22 (A) = 


®22 — A 


det(H 33 - A) 
det yl 3 i + o,i 3 A 


L-(7) 


det(AL 33 — A) 
044 — A 


7 00 (7) 


det(AL 33 — A) 


ci 2 (A) = 
c 2 i(A) = 

C23(A) = 


012 


det(H 33 - A) 
«21 


T°°(7) 


det(H 33 — A) 
det A 32 — 0,23 A 


det(H 33 — A) 


7 / 00 (7) 


7 / 00 (7) 


7°°(7) 


Since | det(AL 33 — A)| > q*|A + 2/0 2 for all A ^ TZ, there exists a positive constant C 4 such 
that 


c ij( A) ^ 


C 4 


| A T k\ T K 
Hence, by virtue of ill. 1(11) . we see that 

1101 II7 2 (7) + 1102117 2 (7) < 


for A ^ 71, i = 1, 2 , 3; j = 1, 2 . 


CV, 


| A + 7Cq| + k 0 


n||7 2 (7) + 11^21| 7 2 (7) + ||s||7 2 (7)) 


whenever Re A + a 0 > 0. Putting these together, we conclude that if Re A > —a 0 then for 
any (r 1 ,r 2 ,s) G (L 2 (/)) 3 equation (14.711 has a unique solution (0!,0 2 ,-0) and it satisfies 
the estimate 


(4.41) 


1101117 2 (7) + || 021| 7 2 (7) + ||0||7 2 (7) < 


c. 


| A + A 0 1 + Nq 


Therefore, we have shown that — C 2 is sectorial in the sector {A G C | |arg(A — k 0 )| < 
7 t/2 }. Once this is established, it is standard that the angle of —£ 2 is less than 7 t/2 (see, 
e.g., pp. 55-56 of [29]). This completes the proof of Lemma [4.21 □ 

4.2. (£o>AL) stability. 

The initial-boundary value problem (12.8(1 is known to have a unique local-in-time classical 
solution (ui(t,x),U 2 (t,x),v(t,x)) if the initial data (w?, u®, v°) is sufficiently regular, e.g., 
u°j G P“(/), v° G C 2 +a (I ) and dv°/dx(0) = dv°/dx(l ) = 0 (see, e.g., [24] [Theorem 1, p. 
111]). We may interpret this as a mild solution of the integral equation on the Hilbert 
space X = (L 2 (/)) 3 as follows: Let <L(t) = r (0i(t, -),0 2 (t, •)). Then 


(4.42) $(t) = e t£2 $(0) + / e (t - T)£2 i/($(s)) dr, <f>(0) G (C(/)) 2 x C\l), 
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(4.43) 


where we define the mapping H : (L°°(/)) 3 —* (L 2 (/)) 3 by 

(ri{<f> i, 02 , 0 )\ 

#($) = 0(01,02,0) , 

\ 001 , 02 , 0 ) ) 

ri(<h) = /i 0 si + 01 , m 2 + 02, fi + 0,0 - aii(x) 0 i - ai 2 (o;)02 - a 13 (x) 0 , 
r 2 ($) = fi (fii + 01 , U 2 + 02, fi + 0,20 - a 2 i(x)0i - a 22 (00 2 - a 23 (x)'ip, 
, s($) = 5 -(mi + 0 i,fi 2 + 02, fi + 0,0 - a 3 i(a;) 0 i - a 32 (x) 0 2 - a 33 (x) 0 . 

By virtue of (14.61) of Lemma 14.21 we have 

t 

||$0)||x < Me _K1< ||<f)(0)||x + [ Me _K1 ^ _T 0|//( < f)(r))||x dr. 


Since |if($)| <C'(|0 1 | + |0 2 | + |0|) 2 , 

(4.44) 

t 

ll*(f)llx < Me—‘I|t(0)||x+CM [ e— “- T )(ll0 1 Mlll. m + II^M||l. m + ||V-(T)|||. (I) )<iT. 


We consider the ODE subsystem, too: 

<90i 

( 4 ‘ 45 ) rS„ 

= a 2 i(x)0i + a 22 (x)cj) 2 + P2(0i, 02,0), 


= an(x)0i + ai 2 (x)02 + Pl(01, 02,0), 


<9f 


where 


Pl (01, 02, 0) = n (0i, 02, 0) + d 13 (x)0, 

P2(01, 02, 0) = r 2 (0i, 02, 0) + 023000, 

Since this is a system of ordinary differential equations with constant coefficients for each 
x, the solution satisfies 


(£(£*)) = exp 033 W) 


01 ( 0,0 
02(0, x) 


exp 


( ( i - T^t*)) (pi{<h(r,x) Ur,x) i,(r,x)^ dT 
K ” 1 P2(01(T,X),02(T,X),0(T,X)) I 


Moreover, by the assumptions on the matrix v4 33 (x), 


exp (L4 33 (x)) 



(0 

R 2 

w 


where || T (pi,P 2 )||m 2 = |pi| + I/O 2 1- Hence, we have a pointwise estimate 

|0i0,0)| + |020,0| < Me _Kf (|0i(O,x)| + |0 2 (O,x)|) 

t 2 

+ M f e~ K{t ~ T) ^2\ P00i(t,x),0 2 (t,x),0(t,x))| G?T. 

n J=1 


(4.46) 
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We now turn to the proof of Theorem 12.41 In order to control the nonlinear terms, we 
first estimate the H l (I) norm of if, which gives us a bound on the L°°(J) norm of if by 
virtue of the Sobolev embedding theorem. 

Lemma 4.5. Let if) = be the solution of the initial-boundary value problem 

for x e I, t > 0, 

for x G dl , t > 0, 
for x G /, 


dib d 2 P 

(4.47) — - D-fop = a 3 i(x)0i + a 32 {x)(p 2 + a 33 (x)ip + s(t,x) 

(448) t^> = ° 

(4.49) ^(0, x) = tpo(x) 


where s 6 CQO,T]; H l (I)) flC 1 ((0,T];L 2 (/)). Then there exist positive constants k and 
Mi such that 


(4.50) 

t 

Wi.-)ll?,.( 7 ) + J e-^dr J 

0 I 




Proof. Let k be a positive number to be determined later. We multiply both sides of 
(14.47P by e kt 'f>(t,x ) and then integrate over (0,£) x I. Observing that fj'ft = 2 _1 ('0 2 )f and 
f’f’xx = (f’'f’ x )x — (f’x) 2 , we obtain after integration by parts: 


n kt 


if(t, x) 2 dx - f (x) 2 dx - [ e kT dr f ip(r,x) 2 dx 


(4.51) + D I e kr dr I \if x (r, x)\ 2 dx 


e kT dr / a 3 3 (i)^(r,a:) 2 + E a 3 j(x)(fj(T, x)^{t, x) + s(r, x)iJj(t, x) j dx. 

3 = 1 


Similarly, we multiply both sides of (14.47)) by e kt i/; t {t, x) and then integrate over (0, t ) x I. 
Noting that = ('4’t'4’x)x — 2 _1 (^) t , we obtain after integration by parts the following: 

(4.52) 

t t 


e kT dr / ^(t, x) 2 dx + 


D 
~2 L 


„kr 


jp x (T,x) 2 dx 


T=t 


- T—0 


Dk 

~Y 


e kT dr / ip x (T, x) 2 dx 


t 



o 


a 33 x) 2 dx + 



a 33 (a;)'?/’(T, x) 2 dx 


T=t 


T —0 



0 
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We use the well-known inequality \ab\ < ~(ea 2 + e l b 2 ) valid for arbitrary £ > 0 
estimate the last two terms of (I4.5ip and (14.521) as follows: 


(4.53) 


e kr dr J ^ a 3j -(a;)0j(r, x) + s(r, x) dx 


o i 3 

t 

— 7 , [ e kr dr f x) 2 dx — R{t ) 


and 


(4.54) 


e kr dr / I a 3 j(x)(pj(T, x) + s(r, x) ) V^(t, x) dx 

r ' 7=1 ' 


0 I 3 

t 

< E — f e kr dr f ijj t (T,x) 2 dx + ^R(t), 


where 


n j J 1 


From (I4.5ip and (14.53)1 follows that 
(4.55) 


e kt / ^j(t, xY dx + D j e kT dr / ^> x (t, x) 2 dx + / e kT dr / |a 33 (a;)|^(r, x) 2 dx 


<— f ^o(x) 2 dx + — f e kr dr j -0(r, x) 2 dx + £ f e kr dr f ^(t,x) 2 dx + —R(t). 


o / 

t 

^kr 


0 


0 


From (I4.52p and il 1.5 Hi it follows that 


e fcT dr f ^ t {T,x) 2 dx + ye fc * f i(j x (t,x) 2 dx + f \a 33 (x)\^{t, x) 2 dx 


(4.56) < y J ip' 0 (x ) 2 dx + i J \a 33 (x)\i/j 0 (x ) 2 dx + ^ J e kr dr J ^ x (t,x ) 2 dx 

II 0 1 

t t 

+ ^ j e kT dr J \a 33 (x)\^(T,x) 2 dx + e' J e kr dr J ijj t {T,x) 2 dx + ^.R(f). 

i oi 


to 


o 
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Adding (14.55)) and (14.561) together results in 


-e kt J x ) 2 dx + — e kt J i^ x (t, x ) 2 dx + (1 — e') J e kr dr J 'i/j t (r,x ) 2 dx 

i i oi 

t t 

+ J e kr dr J i/; x (t, x ) 2 dx + ^1 - 0 J e kr dr J \a 33 (x)\^(r,x ) 2 dx 

0 / 0 I 


n kt 


— ( 7 ^- + £ j f e kr dr j ip(r, x) 2 dx + I |a 33 (a;)|^(t, x) 2 dx 


-\[ [ (x) 2 dx + j |a 33 (a;)|'0o(T) 2 dx + D j ip' 0 (x) 2 dx J + ^ ) R(t). 


Recall from (12.141) that |a 33 (a;)| > 3 k for all x G /. Hence we choose k = £ = min {A, 1} 
to obtain 


1_ l ) / \ a ^ x )\^ T -> x ) 2 dx ~ + £ ) / ^( r , x) 2 dx 


> - ((2 — k)n — k) ip(r, x) 2 dx > 0. 


Also we choose e' — 1/2 and conclude that 


i J ^(t,x) 2 dx + y J ijj x {t,x) 2 dx + ^ J e k{t r) dr J ip t (r,x) 2 dx 

i i oi 

t 

+ j f e - fc (*- T ) dr f x) 2 dx + 0 


< - j / ^o(T) 2 da;+ j \a 33 (x)\ip 0 (x) 2 dx + D j ip' Q (x) 2 dx | e kt 


J i 

t. 


+ {/k +1 ) f e k{t T) dT J ((E a 3j ( x )^( r ’ a: )) 2+s ( T ’ x ) 2 ) dx ’ 

n r i =1 


which implies (14.50)1 . 


□ 
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Put X(t,x) = |0i(f,x)| 2 + \(j) 2 (t,x )\ 2 , Y(t,x) = \^(t, x)\ 2 and Z{t) = \\ 2 Hl{I y So 

far we have established the following inequalities: 

t 

(4.57) X (t , x) <CX M+ cJe-^(X { r, X? + Y { r, x) ) d r, 

0 

(4.58) J(X(t,x) + Y(t,x))dx < C J (X(0, x) + F(0, x)) dxe~ kt 

i i 

t 

+ C [ er k{t - T) dr [ (. X(r,x) 2 + Y{T,x) 2 )dx , 


(4.59) Z(t) < CZ(0)e~ kt + C J e~ k(t ~ T) dr J (X{r, x) + X{t, xf + Y(t, x ) 2 ) dx. 

o i 

Also, by the Sobolev embedding theorem, 


(4.60) 


sup Y(t,x) < CZ(t). 

x£l 


Here and in what follows, we use C to denote various positive constants independent of 
t, X(t), Y(t) and Z(t). From (I4.58|) we get 

t t 

[ e~ Kt - T) dT [ X (r, x) dx < C [ e- k{t - T) e~ kT dr /(X(0, x) + Y(0, x)) dx 


+ C J e~ k(t - T) dr J e~ k{r - a) da j (X(a,x) 2 + Y(a,x) 2 ) dx. 

oo i 

It is straightforward to see that the right-hand side is equal to 

t 

Cfe_ “ / (A ' (o> x)+r(o ' x)) dx+c S ( - t ~ dT I *> 2 + y ^ d * 

I oi 

Since (t — r)e _fc ( t_T ) < Ck'e~ k '^~ T ^ for k! G (0, k), we obtain 

t 

[ e - fc h- T ) dr [ X (r, x) dx 


<C J (X(0,x) + Y(0,x))dx + C J e- kl(t ~ T) dT j (X(r, x) 2 + Y(r, x) 2 ) dx. 

i oi 

Hence, from (14.591) we conclude that 

Z(t) < C(Z( 0) + [ (X(0, x) + y(0, x)) dx) 


+ C g-fc'h-H dr / (X(r,x) 2 + Y(r,x) 2 ) dx 


(4.61) 
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For any measurable set R C /, let 
(4.62) £ R (t)=swpX(t,x). 

xGR 

Then from (14.5711 and (14.6011 . we have 


(4.63) 


Ut) < CUV) + jUt) 2 + c f e- k ^Z(r) dr. 


Decomposing the integral over / into the sum of those over R and / \ R, and denoting 

the Lebesgue measure of R by fi(R), we see 

(4.64) 

/ X(t, x) dx < /i(R)£ii(t) + fi(I \ R)K and / X(t, x) 2 dx </j,(R)£,n(t) 2 + n(I \ R)K 2 , 


where K — U — U_. 

From (14.6lH and '1 1.6 Hi it follows that 

Z(t) < C(Z{ 0) + fi(R)UO) + U \ R)K) 

t 

+ C J e- k ^- T) (fi(R)Ur) 2 + rtl \ R)K 2 + Z(r) 2 ) dr , 

o 

which implies 
(4.65) 

sup Z(t) < C(Z(0)+ M (R)U0)+mV\R)R1+C sup (ix(R)£ r (t) 2 +/li(I\R)K 2 +Z(t) 2 ). 

0 <T<t 0 <T<t 

From (I4.63P we obtain 

Ut) + Z{t ) < c^(0) + cut) 2 + C(Z( 0) + /u(R)UO) + U \ R)K) 

+ C sup {fi(R)Ur) 2 + U \ R)K 2 + Z(r) 2 ), 

0 <T<t 

that is, 

Ut) + Z(t ) < C7(e«(0) + Z( 0)) + C sup (^(r) 2 + Z(r) 2 ) + Cju(I \ R). 

0 <T<t 


Therefore, m(t) = max 0 < T < f (^(r) + Z(r)) satisfies m(t) < C(m(0) + m(t) 2 + fi{I \ R)). 
Now it is easy to show that, for any A satisfying A > max{C', 1}, if we choose £o > 0 so 
small that 

1 + A 2 £q + < —, 

then for any 0 < e < Eo, as long as m(0) < e 2 and / 'j,(I\R ) < £ 4 , the inequality m(t) < Ae 2 
holds for all t > 0. This completes the proof of Theorem 12.41 □ 


Theorem 12.31 on two-component systems is proved in exactly the same way and hence 
we omit it. 

Proof of Corollarv \2.5[ Define a matrix A s (x) by 



/ an (a:) a 12 (x) a 13 (x) \ 

a 2 i(x)/5 a 2 2 (x)/d a 23 (x)/8 
\ « 3 i(^) a 32 (x) a 33 (x) J 
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Then detA 5 (a:) = <5 _1 det A(x), detA 31 (x) = h _1 det Au(x), detA^rc) = detA 22 and 
det A^rc) = h _1 det A 33 (x). Moreover, as 5 j. 0 we have trA s (x) = h _1 a 2 2 + 0(1) and 
tr A 22 (x) = 8 ~ 1 a 2 2 (x) + 0(1). 

Therefore, it is easy to check that under the assumptions of the corollary, conditions 
(I2.1ip - (l2.14p are satisfied for S > 0 sufficiently small. □ 


Proof of Theorem \2. (A To prove Theorem 12.61 we show the conditions of Corollary 12.51 
are satisfied. The only conditions left to check are (I2.18l) - (l2.20p . First we differentiate the 
identity / 2 (wi, u 2 (ui, v), v) = 0 with respect to u\ and v at (ui,v) = (ui,v), respectively, 
and obtain 


(4.66) 


a 21 + a 22 


du* 2 

dui 


0 


du* 2 

and a 23 + cl 22 —— 
ov 


0 . 


Also, by differentiating fi(ui, u 2 (ui, v), v) and g(u\, u 2 (u\, v), v) with respect to u\ and v 
at (u\,v) = (ui,v) respectively, we get 

(4.67) 

rhi* rhi* chi* chi* 

bii=a 11 + a 12 -^, 612 = a 12 ^ + a 13 , 621 = a 31 + a 32 ^, b 22 = a 32 ^ + a 33 . 

OU 1 OV OVb \ ov 


Then we have 


(4.68) 


det A 33 — ana 22 — ai 2 a 2 i 


— CL 22 I On + a 12 


dv^ 

du\ 


= 022^11 > -cia 22 (by (J 22 D)- 


Similarly, by using (I4.66P and (14.671) we see 
(4.69) det An = a 22 &22 > 0 (by (3) and flZZZD)- 

It is straightforward to verify the formula 


(4.70) det A — a 22 (bnb 22 — b\ 2 b 2 i) — a 22 det B, 

which is negative because of (12.71) and (2). □ 


4.3. Quasi-steady state reduction of the ODE subsystem. 

In this subsection we address the question of whether Diffusion Driven Instability can be 
investigated or not based on the quasi-steady state approximation. For 5 > 0, consider a 
system of type 


(4.71) 

(4.72) 

(4.73) 


du[ 


dt 



fliu^U^V 6 ), 

-- f2(u{,U 5 2 ,V S ), 


dv s 

~dt 

dv 5 

dx 


D g^ + 9( u i’ u 2i A 
0 (x <G dl). 


(4.74) 
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Let (ui,U 2 ,v) satisfy 

(4.75) fi(ui,u 2 ,v) = f 2 (ui,u 2 ,v) = g(ui,U 2 ,v) = 0, 

d fo 

(4.76) ^—(ui,U2,v)<0. 
ou 2 

Then by the implicit function theorem there exists a smooth function u 2 (ui,v) defined in 
the neighborhood J\f = {(«i, v) \ \ui — U\\ < rj, \v — u| < rj} of (u\,v) such that 

f 2 (ul,U 2 (ul,v 5 ),v s ) = 0 for (uf, v s ) E J\f and u^ui, v) = u 2 - 


Hence, (ui,u 2 ,v) is a constant steady state of (j4.7ip ~ fj4.73jl . and (ui,v) is a constant 
steady state of its quasi-steady state reduction: 

du\ 


(4.77) 

(4.78) 

(4.79) 


dt 

dv° 

~dt 

dv° 

dx 


= / 1 «^«u°),u°), 
d 2 v 


= D 


dx 2 


+ s(<^K,AA 


= 0 (x E dl). 


Proposition 4.6. Assume that (14.751) and (14.761) are satisfied. If the reduced system 
(I4.77p ~ fj4.79j) exhibits DDI at the constant steady state (tq, v), then there exists a positive 
numberd* such that the unreduced system (I4.7ip ~ fl4.74j) exhibits DDI at (ui,u 2 ,v) as long 
as 0 < 6 < 5*. 


Let B = ipij)i<i,j< 2 denote the Jacobian matrix around the equilibrium (u\,v) of the 
kinetic system for the reduced system, i.e. (I4.77p - fj4.78p with D = 0. Let A 6 = (p4j)i<i,j<3 
denote the Jacobian matrix around {u \, u 2l v) of the kinetic system for the unreduced 
system, i.e., (I4.7ip ~ d4.73l) with D — 0. We prove this proposition by demonstrating the 
following two lemmas, both of which are algebraic in nature. First, in Lemma 14.71 we 
prove that if 5 is sufficiently small then two of the three eigenvalues of A s remain in the 
neighborhood of the two eigenvalues of H, while one of the eigenvalues of A s tends to — oo 
as 5 — * 0, which in particular shows that (mi,m 2 ,u) is stable under spatially homogeneous 
disturbance, provided that 6 is sufficiently small. Second, in Lemma 14.91 we show that 
DDI of the reduced system (I4.77p ~ (l4.79p implies instabilty of (ui,u 2 ,v) for all D > 0. 
Combining these lemma proves Proposition 14.61 

Lemma 4.7. Let B and A s be as above. Assume that df 2 /du 2 (ui,u 2 ,v) < 0. Let Ao,i 
and A 0j 2 be two eigenvalues of B. Then 

( 1 ) if Ao,i 7 ^ Ao ,2 then for each j — 1 , 2 there exists an eigenvalue \$j of A 6 such that 
lim^o | A s ,j - Xo,j\ = 0 for j = 1, 2; 

( 2 ) if Ao.i = Ao ,2 then there exist two eigenvalues A s,j, (j = 1, 2 /, of A 5 such that 
lim^o | A,y - A 0 ,i | = 0 for j = 1 , 2 ; 

(3) there exists an eigenvalue Ay 3 of A 5 such that, Ay 3 is real and lim^o Ay 3 = — oo. 

Proof. Assertions (1) and ( 2 ) are proved by applying Rouche’s theorem. We put aij = ajj 
and A = (dij)i<i,j< 3 - Then afj = for i — 1,3 and j = 1,2,3, whereas af 23 = a 2 j/5 
for j = 1,2,3. Let 4)(A) and T(A) be the characteristic polynomials for B and A s , 
respectively: 

$(A) =A^ — ( 6 n + b 2 2 )X + det B, 

T(A) — — A 3 T (flu T S ^a 2 2 T ® 33 )A^ — (J det An T det A 22 T S 3 det A. 33 )A T S ^ det A. 
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We put 71 = du 2 /dui(ui,v) and 72 = du 2 /dv(u\,v) to see that 

bn = an + 7 i°i 2 , ^12 = 72 °i 2 + ai3, b 2 \ = a^i + 71032, 622 = 72^32 + 033, 

021 + 7 i a 22 = 0 and 72 O 22 + 023 = 0 . 

Using these relations, we obtain (14.67(1 . (14.68(1 and (14.6911 . which yield h'I'(A) = a 2 2 < h(A) — 
t)(A 3 — (an + a 33 )A 2 + det A 22 \). Recall that a 22 < 0 by assumption (14.7611 . Hence, the 
characteristic equation for A 6 reads 

(4.80) a 2 2 < h(A) — 8 if(\) = 0, where if(X) = A 3 — (an + a 33 )A 2 + det H 22 A. 

Let Aon and Aog be the roots of $(A) = 0. Fix a disc D p = {|A| < p} which contains 
both Ao,i and Aog. We apply Rouche’s theorem on D p to show that (14.80(1 has exactly 
two roots Ayi and Xsp inside D p for sufficiently small 5. Moreover, if Ao,i A Aog then, 
by passing to a small disc around A 0 ,, we can prove that Aa ,• —> A 0 ,• as 5 —> 0 for each 
j = 1,2. This verifies (1) and (2). 

In addition, one can easily fold that (14.80(1 has a real root Ay 3 = A _ 1 a 22 + 0{5~ 2 ^) as 
8 —> 0, proving (3). Therefore, the assertions of the lemma hold true. □ 

Remark 4.8. Since det A s = d " 1 det A, if A s is a stable matrix for 8 > 0 sufficiently small, 
then det A < 0 and hence det A s < 0 for all 8 > 0. 


Lemma 4.9. Under assumption (14.7611 . if the reduced system (14.77(1 - ((4.79(1 exhibits DDI 
at (lli,v), i.e., 

(4.81) 6 n >0, ti B < 0 and det£> > 0, 

then (lIi,U 2 ,v) is unstable as a steady state of the unreduced system (14.7111 - ((4. 7411 for all 
D > 0 and 8 > 0. 


Proof. In a way analogous to the proof of (14.68(1 . we obtain 


(4.82) det 4 3 = r 1 ( an ai2 ) = S ~\22611 < 0, 

\a2i a 22 J 

by virtue of our assumptions. Hence H 33 has a positive eigenvalue A* for any 8 > 0. 
Therefore (ui,U 2 ,v) is unstable if D > 0 by the result in §2.1.2 of |9j. □ 


Remark. If we assume, instead of (14. 77D — (14.78(1 . that both species diffuse and consider 
the system 


(«?)* = D !(«?)„ + AK, n u ), u u ), ( V y = D(v") xx + g(ul ^ U ) W U ) 


then the constant steady state (u\,v) is destabilized only for D > D CI with D cr being a 
positive number. This is quite different from the case D 2 = 0 in which (u\,v) is always 
unstable for any D > 0. The same observation applies to the three component system 
(I4.7ip (14.73(1 . See [2] for the case where all three species diffuse. 


5. Proof of statements concerning the example system 
First, we prove uniform boundedness of solutions and classification of different branches. 
Proof of Lemma \3.1\ To obtain bounds on u and v, we see that for v > 0 

(5.1) - (1 + v)u < f r (u, v) < ^-1 + mi U ' 

The hrst inequality of (15.1 p proves that u remains positive if v > 0. 





































32 


STEFFEN HARTING, ANNA MARCINIAK-CZOCHRA, AND IZUMI TAKAGI 


Define 


h(u) := —l + mi¬ 


ll 


1 + ku 2 

If mi > c l\fk then h(u) = 0 has two real roots, and the larger one is u = U 2 + 

\J (mi ) 2 — 4k). Then u(0,x) < U 2 implies u(t,x ) < U 2 for all t > 0, since f r (u,v ) < 0 for 
u > U- 2 - To obtain boundedness of v, note 


(5.2) 


~d 2 u . . ^d 2 v . r^d 2 v m 2 

~ ( " 3 + u)v - + " 3 ” + T 


for u > 0 , proving the second inequality due to the maximum principle. 

If rn 1 < 2 \fk, then max /i(«) < 0. Choose a 7 ^ satisfying 0 < 7 ^ < — ma xh(u). Then 
f r (u, v ) < — 7 /iU for all u > 0 and v > 0. Therefore du/dt < — 7 hU, proving that u(t, x) —$■ 
0 uniformly as t —* 00 . Since g r < — n 3 v + m 2 u 2 / (1 + u 2 ) and m 2 u(t, x) 2 /(l + u(t, x) 2 ) —> 0 
as t —> 00 , we conclude that lim^oo v(t, x) = 0 , uniformly on /. 


□ 


Proof of Lemma EJI To obtain uniform boundedness of solutions, note 


du[ 

dt 

.du 2 

’if 

dv 5 
~dt 


> {61 


u 1 


1 1 + k(u{) 2 

> —(^2 + a ) M 2 ; 

d 2 v 5 


V\ 


pv s )u{, 


> 7' 


<9x 2 


(1/3 + / 3 uJ)i; 4 , 


<9(m 3 + hu: 




= — ( z 7?4 + ^2^2) + $1 


(u 


S\2 


(5.3) 

(5.4) 

(5.5) 

and for <5 < 1 

(5.6) 

Consequently, 

(5-7) .■/ x- 

t-y 00 mm(/ 7 , v 2 )k 

In particular, if wf(0, x) + Su 2 (0,x) < #i/(Acmin(i 7 , v 2 )), then uf(t,x) + Su 2 (t,x) < 
#i/(fi;min(z 7 , u 2 )) for all t > 0. Therefore, if u 2 ( 0, x) < 6i/(5nmm(vi,v 2 )) then we have 


1 + k(u{) 2 


< - min(z/i, v 2 )(u[ + Su 2 ) + 


0 1 


K 


lim sup (u\ + Su s 2 ) < 


0 i 


(5.8) 


" s ,if 




+ 


09 


K 


k 8 min(z 7 , z/ 2 ) 

yielding the result for v s . 

Now assume Q\ < 2uiy/ti and put h\(£) = —Vi + 6 h£/(l + k£ 2 ). Then max hi (£) 
is negative, and hence we can choose 71 satisfying max hi (£) < —71 < 0. Therefore, 
—Viu\ + 9i(u{) 2 /(l + k{u\) 2 ) < —71 u{ for u\ > 0. We thus hnd from (j5.6[) that d{u\ + 
Su 2 )/dt < —('fiu{ + v 2 u 2 ), which yields the uniform convergence uf(t,x) + Su 2 (t,x ) —» 0 
as t —* 00 . Once this is established, we have dv 5 /dt < n fd 2 v s /dx 2 — u 3 v s + c(f) with 
Hindoo c(f) = 0 as in the proof of Lemma [3. 11 yielding the result for v s . □ 


Proof of Lemma HOI If (u,v) is a steady state of system (j3.3h - fj3.5jl . it satishes 

du zi? 

(5.9) — = f r (u,v) = -(l + r)« + mi i + fa2 = 0. 
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This is satisfied if and only if either u = 0 or 

'll 

(5- 10 ) -(l + ^+rrat = 0. 

Solving (I5.10P for u yields 

(5.11) u±(v) = ^YTv) ( mi ± ]/ m i~ 4k ( 1 + v ) 2 ) 

proving Lemma [3.31 □ 

The spatially homogeneous steady states are defined as common roots of two polyno¬ 
mials (1 + ku 2 )f r (u,v) = 0 = (1 + ku 2 )g r (u,v). We therefore want to distinguish the 
different branches by a criterion which is easier to handle and note that the branches of 
steady states can be characterized in larger generality: 


Lemma 5.1. The solution u(v) of f r (u,v) = 0 has three branches: 

(5.12) u 0 (v) = 0, 

(5.13) ' ’ 1 


“±M = 2 k(l+v) ( m ‘ ± \A n ?- 4fc ( 1 + ' < ’) 2 ) • 


2 k(l + v) 

For u 7 ^ 0 and 0 < v < v r — m\/ (2 \fk) — 1, the following are true: 

( 1 ) u{v) = u + (v) if and only if < 0 ; 

(2) u{v) = U-{v) if and only if > 0. 

Proof. Note that 

1 


(5.14) 


u + (v) = 


2 k(l + v) 


m i + \J m\ — Ak(l + v ) 2 


implies ^ u + < 0. Combining this with 


(5.15) 
yields 

(5.16) 


n_i_u_ = 


k' 


d_ 

dv 


u _ > 0 . 


□ 

Using this characterization, we can prove existence of spatially homogeneous steady 
states and classify them: 

Proof of Lernrna 1 3 . 51 Clearly, ( uq,vo ) = (0,0) is a spatially homogeneous steady state. 

To obtain nontrivial spatially homogeneous steady states, observe that a solution (u, v ) 
is a spatially homogeneous steady state only if 


(5.17) 

Solving f r (u,v ) = 0 yields either u = 0 or 


= /„(„,„) = -(l +v) u + miT f 1 - 2 


= 0 . 


(5.18) 

Inserting (I5.18p into 

(5.19) 


u 


= — (1 + v). 


1 + ku 2 mi 


dv u ^ 

- = 9r (u,v) = + u)v + m 2T - il - 2 


= 0 
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and solving for v yield 

(5.20) v Ir ,, r (u) 

Note that Vf r:9r (u) is A:-independent. 


m 2 u 

mi h 3 + (1 - ^)u 



Figure 5.1. Nullclines of f r and g r and Vf r)9r for parameters mi = 
1.44, m 2 = 2, n 3 = 4.2, k = 0.1. 


The nullclines of f r and g r for u ^ 0 are described by 

'll 

(5.2!) v fr=0 (u) := -1 + mi^, 


(5.22) 


Vgr=oi u ) 


m 2 u 2 
/13 + MI + ku 2 


If there exists u > 0 such that v 9r= o(u) = Vf r , 9r (u), then {u,Vf r , 9r (u)) is a homogeneous 
steady state. 

Note that Vf r=0 (u ) has a unique positive maximum at u = 1/Vk, is strictly concave on 
[ 0 , y/3/y/k), strictly convex on (y/3/Vk, 00 ) and satishes 


(5.23) limu /r=0 («) = lim v fr=0 (u) = -1. 

u \U u/'oo 

Recall rrii < m 2 . Defining l{u) := m 2 w 2 /(/i 3 + u), we see that 

(5 ' 24) ' (u) " v *-° (u) = T+W) 

(5.25) Z(u) - v 9r=0 (u) > 0, 

(5.26) lim ||l - v 9 r= 0 || L oo([o,c]) = 0, 


u 


5 


on any finite interval [0, c]. Note also that Vf r , 9r {u) is strictly increasing and continuous on 
M>o \ {(mi/is)/(m 2 — mi)}, nonnegative and convex on [0, (mi/i 3 )/(m 2 — mi)), negative 
on ((mins) /(m 2 — mi), 00 ) and satishes 


(5.27) 


lim 1 (u) = 00 , 


11 / 

x m 2 —m^ 


lim Vf r<9r (u) = - 00 . 

m l M3 •> 1 1 ' / 


(5.28) 
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Furthermore, it holds that 

(5.29) 

(5.30) 


v fr,gM = %.=o(0) = 0, 


dv Br= 0 


( 0 ) = 0 < 


dv 


f 

JriQr 


( 0 ). 


du du 

Hence, u 9r . =0 (c) < Vf r: g r (e) for £ > 0 small. 

Combining this with (15.27(1 . the uniform boundedness of v 9r= o and the strict convexity 
of Vf r ,g r on (0, (mifis)/(m 2 — mi)), we see that if ( u , l(u)) and uvfr,g r ( u )) intersect twice in 
the region 93 = (0, (mi/i 3 )/(m 2 —mi)) x (0, oo), then there exists k\ such that (u, v 9r=0 (u)) 
and (u, Vf r ^ r (u)) intersect twice in the same region 93 for all k < k\. Moreover, if ( u , l(u)) 
and (u,Vf rt g r (u)) do not intersect twice, neither do (u,v gr=0 (u)) and (u,Vf rt g r (u)) for any 
k > 0. Now we solve the ^-independent equation l(u) = Vf rt9r (u), i.e., 

9 

m^i 1 

(5.31) 


7T12U 


m 2 


u 


h3 + U mi J u 3 +(l- =*) 


u 


for u and obtain for 

mi 


(5 ' 32) U± 2(m 2 -m 1 ) 

Clearly, u± > 0 holds if and only if 

(5.33) 


± 


miii 3 


h 3 


2(m2 — mi) J m 2 — m\ 


(5.34) . . 

\ 2 (m 2 — mi) 

Observe that (I5.34p is equivalent to 


miii 3 


A 6 > -, 

mi 

2 


> 


h3 


m 2 — mi 


(5.35) 

Recall m 2 > mi and note that 


9 mi — 2m.9 1 

fl 3 + 2-5- Us -f- 2 > 0’ 

"— 1 mf 


9 

mf 


m 2 

mi 


m 2 

mi 


> 0 


Hence, (15.35(1 is satisfied if and only if 


(5.36) 


h 3 < — 

mi 


2 m .2 — mi 
mi 



— 2 a / ( —- 1 - 1 


or 


(5.37) 


/ 1 3 > 


mi 



Inequality (I5.36jl is never satished if ((5.33(1 holds, because it holds that 

2 m -2 — mi 


(5.38) 


1 

mi 


m 2 


mi 


m i 


m 2 


- 2 t — <—, 


mi 


mi 


implying /x 3 < 1/mi if (15.36(1 holds, which contradicts (15.33(1 . 
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On the other hand, since m 2 > mi, it holds that 

2m2 — mi 


(5.39) 

Therefore, if 

(5.40) 

(5.41) 


+ 2 , 


m\ 


h 3 > — 

mi 

m2 > mi, 


2m2 — mi 
mi 


m 2 

mi 


+ 2i 


m 2 

mi 


> 1 . 


m 2 

mi 


both (5.33) and (5.34) are satisfied, and hence (■ u,l(u )) and (u,Vf t9 (u)) intersect twice in 
94. This establishes the existence of two homogeneous steady states for sufficiently small 
k. 

To prove that both steady states are of type w_, we use the fact that the steady states 
component u is less than (mifi 3 ) / (m 2 — mi). Since Vf r=0 (u) = —1 + miu/(l + ku 2 ) is 
strictly increasing on (0, 1/y/k), it holds that at any steady state u* 


(5.42) 


dv(u) 

du 


(u*) > 0 


for \fk. < (m 2 — mi)/(mi/], 3 ). By the characterization Lemma f5.ll both spatially homo¬ 
geneous steady states are of type (u-(v),v) for k < ((m 2 — mi)/(mifi 3 )) 2 . □ 


Existence of spatially inhomogeneous steady states can be proved similarly to [2U] or 
[TH] . To apply this work, we solve f r (u,v) — 0 in u and substitute this into 

u 2 

(5.43) g r (u,v) = -(fjL 3 + u)v + m 2 l ku2 . 

The shape of g r (u(v),v) for different branches of the solution u(v) of f r (u, v) = 0 is shown 
in Fig. 15.21 



Figure 5.2. Illustration of the right-hand side of —d 2 v/dx 2 = g r (u,v) for 
different branches of the solution u(v) of u t = f r (u, v) = 0. The parameters 
for illustration are D = l,mi = 1.44, m 2 — 2, fi 3 — 4.2. We can observe 
that all nontrivial homogeneous steady states are of type w_. 
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Lemma 5.2. Assume 2 \fk < m\ < m 2 . Then the following inequalities hold: 

(5.44) <7 r (u±(0),0) >0, 

(5.45) g r (u + (v), v) > g r (u + (v r ),v r ) = g r (u-(v r ), v r ) for 0 < v < v r , 

(5.46) g r (u + (v), v) < g r (u + ( 0), 0) for 0 < v < v r . 

If there exist two positive spatially homogeneous steady states of type U-, it holds that 

(5.47) g r (u±(v r ), v r ) > 0, 

(5.48) g r (u + (0), 0) < 00 , 
and 

(5.49) ^g r (u-(v),v) -> — - /x 3 + 2 ^ 2 2 — (1 + v), 

dv nii m \ 

uniformly on any finite interval v G (0, c) as k \ 0. 


Proof. Inserting v = 0 into the right-hand side (5.43) leads to 
(5.5°) Sr MO), 0 ) = m 2 1+ M ^ (q)2 > 0. 

yielding (15.440 . 

We differentiate g r (u + (v),v) with respect to v and obtain 


(5.51) 


d_ 

dv 


■9r(u+(v), 


V ) 


(1 + v)m\ + 2k {1 + n) 3 mi/i 3 
2 h(l + v) 3 mi 


< 0 


mf + 4 fc(l + v) 3 (777,2 — 777 i) 


for v > 0, hence (15.45[) and (I5.46P hold true. 

If there exist two spatially homogeneous steady states of type u_, then g r (u±(v r ), v r ) 
must be positive because g r (u_(0), 0) > 0 and g r (u_(v),v) has exactly two roots of or¬ 
der one in the interval 0 < v < v r . Hence (I5.47P is proved. Inequality (I5.48p follows 
immediately from n + (0) <m\/k and (15.501) . 

To see (15.491) . we differentiate g r (u-(v),v ) with respect to v: 

d 777 3 + 4h(l + v) 3 (?77 2 — 777 1 ) —?77 2 — 2h(l + 7>) 2 ?77i/i 3 

dv 9r u - v ,v mi 2k(l + v) 2 -Ak(l + v) 2 + m\ + 2fc(l + v) 2 mi 

m\ — ?77 2 \J —4h(l + v) 2 + m\ 2(1 + v)(?77 2 — mi) 

— —/i 3 - , = — - = . 

27?7ifc(l + n) 2 y —4fc(l + v) 2 + m\ m\\J —4h(l + v ) 2 + 777 ^ 

The limits for k —> 0 of the first and the last term are clear. To obtain the limit of the 
second term, we use l’Hopital’s rule and obtain 


rK ^ jk( m l - 777! \J 4fc(l + V) 2 + 777 2 ) _ m X 

' ^(2fc(l + u) 2 a/— 4/c(l + V ) 2 + 777?) -6fc(l + u) 2 + ?«?’ 

leading to 

(5.53) lim i r (u_(t;),t;) = — - g 3 + 2 ( m2 ~ mi ) (i + v ) 

V ' k\0dl> V W ' 777i m? 

which proves (15.491) . 


□ 
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Using this knowledge, especially about the sign of g r , we can prove existence of steady 
states with jump discontinuity. The proof is oriented on the shooting method and follows 
the principle presented in ra, a and [18] . However, for completeness we state 

Proof of Lemma fXhl The concept of the proof is illustrated in Figure 15751 Step 1. We 
construct steady-states of (I3.3l) - (l3.5p by using two building blocks. The first one is the 
solution of the initial value problem 

v"(x) = —g r (0,v(x)) = gsv(x) for x > 0 , 

(5.54) t/(0) = 0, 

r>(0) = a, 

which is solved explicitly: 

(5.55) v(x) = a cosh( A //X 3 x). 

We denote this solution by v(x\a). The second one is the solution of the initial value 
problem 

w"(x) = —g r (u + (w(x)),w(x)) for x > 0 , 

(5.56) w'(0) = 0, 

w( 0 ) = b. 

We denote this solution by w(x,b), where the initial data Wq is taken from the interval 
(0, v r ). Denote — g r (u + (w), w) by h(w), which is dehned for w E [0, v r ]. Then we know 
from Lemma 15.21 that 

(5.57) h( 0) < h(v r ) < 0, 

(5.58) ~r~( w ) >0 for w E [0, ry). 

aw 

Clearly w(x) is positive and monotone decreasing in a certain interval 0 < x < x b where 
w(x b ) = 0. Note that x b < y/2b/\h(b)\ since w(x) = b + h(w(6x))x 2 /2 with 0 < 9 < 1. 

Step 2. We would like to connect w(x]b) with v(x-,a ), or v(x',a ) with w(x]b), so 
that the boundary conditions are satisfied. Fix a number c in the interval (0, b) and let 
x c = x c (c, b) E (0, Xb) be such that w(x c ] b) = c. Now we connect w(-; b) to v(-; a) at w = c 
as follows: We would like to choose an appropriate a E (0, c) and a y c — y c {c, a) > 0 so 
that the matching conditions 

(5.59) v(—y c \ a) — c and v'(-y c ; a) = w'(x c \ b) 
are satisfied. This is possible if and only if 

(5.60) \w'(x c (c,by,b)\ < 

is satisfied. Indeed, if (I5.60P holds, then there is a unique £ c < 0 such that tanli£ c = 
w'{x c \b) / (y/jlzc). Then we obtain y c = — id and a = c/ cosh £ c . We say that w(x]b) 
is switchable to v(x) at w = c if 05.59P is satisfied. We remark that for each b there 
exists a unique c* = c*(b ) E (0, b) such that w(x] b) is switchable to v{x) at c if and only 
if c*(b) < c < b. This is proved easily by using the following three facts: (■ w'/w )' = 
h(w)/w — (w'/w) 2 < 0 , w'(0)/w(0) = 0 and lim^ Xtuo w'(x)/w(x) = —oo. 

Next, fix a number 7 E (a, v r ) and let y c = y c ( 7 , a) > 0 be such that v(y c ] a) = 7 . We 
extend w(x] b) to — x b < x < 0 by defining w(x] b) = w(—x\ b) for x E [—x b , 0]. This time 
we would like to find b E (a, v r ) and x c = x c ( 7 , b) E (0, x b ) such that 

(5.61) w(—x c ]b) = 7 and w'(—x c ;b)—v'(y c ',a), 
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which is equivalent to the following conditions: 

(5.62) w(x c -, 6 ) = 7 and w'(x c ;b) = -v'(y c ; 7 ). 

We say that v(x; a) is switchable to w(x) at v = 7 if there exist b and x c ( 7 , b) satisfying the 
matching condition (15.611) . Likewise for w(x; b ), we shall prove that there exists a unique 
c* = c*(a) such that v(x; a) is switchable to w(x) at v — 7 if and only if a < 7 < c*(a). 

Step 3 . Assuming the existence of c*(a), we describe how to construct a steady state 
of (I3.3[) - (I3.5D . First we choose b G (0,ry) arbitrarily and then choose any c G (c*(b),b). 
We now obtain a and ?/ c (c, a) which satisfy (15.591) . Let us put 


W(x) 


w(x; b) for x G [0, x c (c, b )), 

v(x - x c (c, b) - y c (c,a ); a) for x G [x c (c, 6), x c (c, 6) + y c (c,a)}. 


Clearly W\(x) = W((x c (c, b) + y c (c, a))x) is a weak solution of (I3.3f) - (l3.5f) for D = l/(x c + 
y c ) 2 . Next, choose 7 G (a, c*(a)) arbitrarily and find a pair (6 l5 x c ( 7 , bi)) satisfying (15.62p . 
We extend W (x) to the interval [0, x c (c, b ) + y c (c, a) + y c ( 7 , a) + y c ( 7 , 61 )] by putting 


f W(x) for x G [0, Li], 

W(x) = < v(x — L\\bi) for x G [Li, L 2 \, 

[u;(x - L 3 ;60 for x G [L 2 ,L 3 \, 

where Li = x c (c, b) + y c (c, a), L 2 = L x + y c ( 7 , a) and L 3 = L 2 + x c ( 7 , 61 ). Actually, IF(x) 
is extended to [0, L 3 + x^J. If we put kF 2 (x) = W(L 3 x), then this gives rise to a weak 
solution of (13. 3 p (13.51) for D = I/L 3 , which has two switching points. We can repeat these 
procedures to obtain steady states of (13.31) — (I3.5D with an arbitrary number of switching 
points. In the same manner, starting with u(x;a), one obtains another type of steady 
states with any number of switching points. 

Step 4 ■ Now we prove the existence of ( 6 , x c ( 7 , 6 )) satisfying (15.62)) . Since w(x]b) is 
not an elementary function, we investigate its qualitative properties. First, let 0 < b\ < 
b 2 < v r . Then w(x m , bi) < w(x\ b 2 ) in the common interval of existence. For, 0(x) = 
w(x; b 2 ) — w(x; 61 ) satisfies 4>" = h'(w(x ; b \) + 6*0(x))0, 0'(O) = 0 and 0(0) = b 2 — b x > 0 
with 0 < 9 < 1 . Since h'(w) > 0, this implies 0(x) > 0 for x > 0 as long as it exists. Next, 
fix 7 G (0, v r ) and let y c = y c ( 7 , a) > 0 be such that v{y c \ a) = 7 . Then for each b G ( 7 , v r ) 
there exists a unique x (7 ; 5) > 0 such that tc(x (7 ; 6 ); 6 ) = 7 . Notice that x( 7 , b) is a 
continuously differentiable function of b in ( 7 , v r ) due to the implicit function theorem, 
since w(x( 7 , 6 ); b) =7 and w'(x{ 7 , 6 ); 6 ) < 0. Moreover, we claim that if 7 < 61 < b 2 < v r 
then w/(x( 7 , &i); &i) > w'(x{p), b 2 ); b 2 ). To verify this assertion, we define 


(5.63) 


W 



0 


Then from ()5.56j) it follows that 

(5.64) (u/(x; b)) 2 = 2 (. H(w(x ; b )) — H(b )), 


and hence 


(5.65) K(x( 7 ; 63 ); & 2 )) 2 - K(x( 7 ; 60; 60) 2 = -2 (H(b 2 ) - H (60). 

Recalling (5.57) and b 2 > 61 , we find that (w/(x( 7 ; 60; 60) 2 > (tc'(x( 7 ; 61 ); 60) 2 . Since 
w'(x; bj) < 0 for x G (0, x&.), j = 1, 2, we conclude that ta / (x( 7 ; 6 2 ); 6 2 ) < w/(x( 7 ; 61 ); 61 ). 
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We have proved that b n- w'(x( 7 , b);b) is a strictly decreasing continuous function. It is 
easy to check that (i) w'(x( 7 , 6 ); b) —» 0 as b l 0, (ii) w'(x( 7 , fe); b) 7 ) — H{v r )) 

as b 'l' v r , and (iii) if we put p*( 7 ) = ^ 2 (^( 7 ) - #(dv)) then ^( 7 ) = h(i)/p*{i) < 0, 
P*(0) = \J-2 H(v r ) > 0, and p*{v r ) = 0. Therefore, given g G (0,p*(7)), there exists a 
unique b G ( 7 , v r ) such that w'(x(^, b), b) = —q. 

Hence, our problem reduces to asking whether \v'(y c (j, a); a)| G (0,p*(7)) or not, 
given a and 7 . Since v(x;a) = a cost^yZ/Lyc), we see by elementary computations that 
xfipiydj, a) = cosh -1 ( 7 /a) and 

v'illci'y, cl)] a) = - 1 // 7 iasinlW 1 (coslW 1 ( 7 /a)) = ~ a 2 . 

These facts imply that there exists a unique c*(a) G ( a,v r ) such that a/A t 3 (c*(a ) 2 — a 2 ) = 
p*(c*(a)) and |n'( 7 / c ( 7 , a); a)| < p*(q) if a < 7 < c*(a), while ^'(^( 7 , a); a)| > ^*( 7 ) if 
a G (c*(a), ty). This proves the existence of c*(a). 

Step 5. Finally we examine the range of D which has monotone decreasing solutions. 
Note that c* is characterized by the property that \w'(x c (c*,b)]b)\ = . Therefore, 

tanh£ c —>• 1, i.e., £ c —> —00 as c | c*( 6 ), so that y c —> +cx) and a —> 0. This means 

that D = l/{x c + y c ) 2 —* 0 as c i c*( 6 ). On the other hand, as c f we have (i) 

x c —y 0 and (ii) w'(x c ; 6) —» 0 and hence —> 0. This implies that y c —> 0 and a —> c 

as c t 0. I 11 particular, D = l/{x c + y c ) 2 —>• +cxd as c "\ b. Observing that x c and y c 

are continuous function of c, we conclude that for each D > 0 and for each b G (0, v r ) 

there exists a monotone decreasing solution Hq (x:b) of (13.3D - (13.51) as stated in Step 3. 
Clearly W\{x\b\) 7 ^ W\(x',b 2 ) if 0 7 ^ b 2l we hence obtain infinitely many steady states 
for each D > 0. (The arguments above also prove the existence of infinitely many steady 
states which are not monotone in x for each D > 0.) This completes the proof of Lemma 

[3761 □ 



Figure 5.3. Illustration of the construction of weak steady states in the 
proof of Lemma 13.61 


Now that we have established the existence of spatially homogeneous steady states and 
spatially inhomogeneous steady states, we turn to the proof of the lemmas stating that 
the conditions of Theorem 12.31 are satisfied: 
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Proof of Lemma \ 3. 7[ First, we calculate the Jacobian matrix of the kinetic system for 
given arbitrary ( u,v ): 


(5.66) 


B t x \ . = An b l2 \ _ (~{l + v) +m 1(1+fc “ 2)2 
\ &21 ^22 . 


-W 


- V + m 2 (1^2)2 -(/i 3 + M)i‘ 


The signs of bij for u = u 0 = 0 follow immediately by inserting u = 0; also the signs of 
b \2 and b 22 are obvious for arbitrary u. 

For f r (u, v) — — (1 + v)u + rriiu 2 / (1 + ku 2 ) it follows 

, ^ . 2m,iu 

b u = -(l+v) + 


(5.67) 


(1 + ku 2 ) 2 

= ~( 1 + ^)+ 1 2 ,M + v ) 


1 + ku 2 


1 — ku 2 


(1 + u). 


1 + ku 2 

Combining (15.6711 with U-(v) < 1/y/k < u+(v) yields the result for bn. 
We investigate b 2 n 

2 u 


b 2i = -v + m 2 


(1 + ku 2 ) 2 


(5.68) 


m 2 . . m 2 

= -v + —{l + v) + — 

mi mi 


-1 — v + mi 


2 u 


(1 + ku 2 ) 2 


m 2 


mi 


m 2 


= -h-1 iH- -bn- 


mi 


m 2] 


mi 


The continuity of b 2 i and &n(ry) = 0 imply the result for b 2 \. 
In the proof of Lemma 15.21 we saw in (15.511) that 

(5.69) ^-g r {u+{v),v) < 0, 


dv" 


and in Lemma 15.11 that 
(5.70) 


A +W < °' 


It is left to prove the dependence of the determinant of B(x) on dg(uf r= o(v),v)/dv. To 
see this, recall from the identity f r (u(v),v) = 0 that 


(5.71) 


, du , „ . d .... du. , 

Oii 3 - + Oi 2 = 0, and — g r (u(v), v) = — b 2 i + b 22 - 
dv dv dv 


The first equation yields du/dv = —bi 2 /bn- Hence, from the second equation we obtain 

d hi 2 , det B 

—g r {u{v),v) = 6 2 i + 622 = -• 

dv b\\ b\ 1 

From the first equation of (15.71[) we get the expression bn = — bi 2 /(du/dv). This in turn 
give us 

det B = bu—g r (u(y),v) 

(du\ 1 d 


= -M*d Tv 9Mv) ' v) ’ 


as desired. 


□ 
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Proof of Lemma [?P1 First, we prove that spatially homogeneous steady states of type 
(u-(v),v) are unstable as steady states of system fid.3D - (13.51) for any D > 0. Lemma 13771 
states 


(5.7 2) 9 u fr |(u_(i;),u) 0) 

proving instability since autocatalysis implies instability, see [l 6 lj. 

Letting B again denote the Jacobian matrix of the kinetic system at the steady state, 
we recall that a steady state of the kinetic system is stable if and only if tr (B) < 0 and 
det B > 0. Hence, we calculate the trace: 

trB= (- 1 + TTk ^} {1 + v) - {ll3 + u) 

2 m 2 u 

1 + ku 2 m i v 

since u(/i 3 + u) = m 2 u 2 / ( 1 + ku 2 ) = m 2 u(l + v)/m,\ due to (5.18) and (5.19). From 

u 

(5.74) v — — 1 + m\ --—- < — 1 + mill, 

1 + ku 1 

it follows that 

(5.75) bn + b 22 < 1 + 2--—--^ (1 + u). 

\ mi —1 + niiii ) 




Therefore, it is easy to conclude that tr B < 0 if m\ < ^/m 2 . 

To prove detH > 0, we recall that g r (n_(0), 0) > 0, g r (u_(v r ),v r ) > 0 and g r (u_(v),v) 
has exactly two roots V- < v + G (0,u r ), see Lemma 15.21 It follows from the proof of 
Lemma 15.21 that 


(5.76) 


dg r (u_(v),v) 

dv 

dg r (u-(v),v) 

dv 


(n+) > d(k), 
(v-) < -c 2 (k). 


Note 


(1) U-(y) > 0, see Lemma [5711 

(2) U-(v) > 1/mi, see proof of Lemma 15.11 

(3) det B = uj^g r {u-(v),v)/^-{v), see Lemma O 

Assuming dg r (u_(v),v)/dv(v±) —> ±|c±| as k \ 0, (I5.76|) therefore implies 


(5.77) 


det B(u-(v + ), n + ) > 0, 
det B{u-{v-),v-) < 0. 


It is left to prove dg r (u-(v),v)/dv(v±) 0 as k —* 0. To see this, recall that (u-(v-),v-) 

and (u_(n + ), v + ) converge towards two different, finite limits as k \ 0 , see the proof 
of Lemma 13.51 Since (15.49)1 states that dg r (u-(v),v)/dv converges uniformly towards an 
affine-linear function on any finite interval, we obtain the bound on the derivative. □ 
Finally we come to the 

Proof of Lemma \3. 1 1\ Assertion (1) follows directly from Lemma [33] since / 2 (rtf, u 5 2 , v s ) = 
0 is uniquely solvable and the rescaling is linear with positive constants. 

To investigate stability of spatially homogeneous steady states, we notice 

(5.78) d 2 f 2 (u[,u 2 ,v 5 ) = ~(n 3 + S) < 0 . 
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First we verify the stability or instability for the kinetic system in assertions (2), (3). Due 
to (15.781) and Lemma 14.71 for all X s F a(A s ) it holds that either X s G I h c ^ B c(S)( X) 
with lim^o c(S) = 0 or X s E M — > — oo, where B c ^)( A) denotes the open disk {/i e C | 
\jji — A| < c(<5)}. By Lemma 13.91 the Jacobi matrix B for (w_(u_), u_) has a positive eigen¬ 
value, whereas both eigenvalues of B for (m_(w + ), ry_) have negative real part. Therefore 
assertions for the kinetic system in (2) and (3) hold true. Second, we obtain the destabi¬ 
lization of spatially homogeneous steady state (m_(u + ), v + ) for D > 0 by applying Lemma 
4.9, completing the proof of assertion (3). Since (w_(n_),n_) is an unstable equilibrium 
of the kinetic system, it is unstable also as a steady state of (13. ip . To obtain stability of 
steady states with jump discontinuity and also of (0,0,0), we show that the conditions 
of Theorem 12.61 are satisfied. Condition (1) is satisfied, see Lemma [3.71 and Corollary 13.81 
Condition (2) is satisfied due to Lemma 13.21 and condition (3) is satisfied due to (15. 78[) 
and 

(5.79) a 33 = -v 3 - Pu{ < 0. 

Assertion (5) is a consequence of Corollary 13.81 and Theorem 12.61 Assertion (4) follows 
from Lemma 13.71 and Theorem 12.61 This finishes the proof. □ 
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